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Cj^ ' Abstract. Recent results strongly challenge the canonical picture of massive star winds: various evidence indicates that cur- 

rently accepted mass-loss rates, M, may need to be revised downwards, by factors extending to one magnitude or even more. 
This is because the most commonly used mass-loss diagnostics are affected by "clumping" (small-scale density inhomo- 
geneities), influencing our interpretation of observed spectra and fluxes. 

^^ , Such downward revisions would have dramatic consequences for the evolution of, and feedback from, massive stars, and thus 

robust determinations of the clumping properties and mass-loss rates are urgently needed. We present a first attempt concerning 
this objective, by means of constraining the radial stratification of the so-called clumping factor 

To this end, we have analyzed a sample of 19 Galactic 0-type supergiants/giants, by combining our own and archival data for 
O ' Hq,, IR, mm and radio fluxes, and using approximate methods, calibrated to more sophisticated models. Clumping has been 

}—{ ' included into our analysis in the "conventional" way, by assuming the inter-clump matter to be void. Because (almost) all our 

C/3 . diagnostics depends on the square of density, we cannot derive absolute clumping factors, but only factors normalized to a 

certain minimum. 

This minimum was usually found to be located in the outermost, radio-emitting region, i.e., the radio mass-loss rates are the 
lowest ones, compared to M derived from H^ and the IR. The radio rates agree well with those predicted by theory, but are only 

^^ ' upper limits, due to unknown clumping in the outer wind. H„ turned out to be a useful tool to derive the clumping properties 

inside r < 3. . SRj,. Our most important result concerns a (physical) difference between denser and thinner winds: for denser 
winds, the innermost region is more strongly clumped than the outermost one (with a normalized clumping factor of 4. 1 ± 1 .4), 
whereas thinner winds have similar clumping properties in the inner and outer regions. 

Our findings are compared with theoretical predictions, and the implications are discussed in detail, by assuming different 
scenarios regarding the still unknown clumping properties of the outer wind. 

Key words. Infrared: stars - radio continuum: stars - stars: early-type - stars: winds, outflows - stars: mass-loss 



1 . Introduction l2005l) . Unfortunately, however, our knowledge of these objects 

is not as complete as we would like it to be, and present ef- 
In the last few years, massive stars (Mzams ^ 10 Mq) have f^^^g concentrate on modeling various dynamical processes in 
(re-)gained considerable interest among the astrophysical com- j^e stellar interior, as well as in the stellar atmosphere (mass 
munity, in particular because of their role in the develop- loss, rotation, magnetic fields, convection, and pulsation), 
ment of the early Universe (e.g. , its che mical evolution and 
re-ionization: lBromm et alJ200ll but also lMatteucci & Calural 

Most important in this regard is the mass loss that occurs 

Send ojfprint requests to: J. Puis through supersonic winds, which modifies evolutionary time- 

* based in part on observations obtained with the VLA operated by scales, chemical profiles, surface abundances and luminosities. 

the National Radio Astronomy Observatory (NRAO) As shown by numerous stellar-evolution calculations, changing 
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the mass-loss rates of massive stars by ev en a factor of t wo has 
a dramatic effect on their evolution ( Mevn et et alJl994l) . 

The winds from massive stars in their 0-, B- and A- 
supergi ant phase are well d escribed by radiation-d riven wind 
theory JCastor et all 1 19751 Pauldrac het al.lll986l) : the even 
stronger mass outflows observed during their Wolf-Rayet 
(WR) and Luminous Blue Variable (LBV) phases are also 
thou ght to be driven by radiation pressure (for recent progress, 
see Grafener & Hamann 2005 for the case of WRs and 
lOwocki et al. 2005, for LB Vs). 

Notwithstandi ng its considerable succe sses (e.g., 



iVink et a l. 2000; Kudritzki 2002; Puis et al. 2003), the flieory 
is certainly over-si mplified. Stellar rotation (e.g.. .Owocki et alJ 
I1996HPu1s et alJl9 99 and references therein), and the intrinsic 
instability of the line-driving mechanism (see below), produce 
non-spherical and inhomogeneous structure, observationally 
evident from, e.g.. X-ray emission and lin e-pr ofile variability 
(for s ummaries, see lKudritzki & Pulsl2000l and lOskinova et alJ 
I2OO4I regarding the present status of X-ray line emission). 
As long as the time-dependent structuring of stellar winds 
is not well understood, we cannot be sure about even their 
"average" properties, such as mass-loss rates and emergent 
ionizing fluxes. Even worse, most spectroscopic analyses of 
hot stars aiming at deriving stellar and wind parameters have 
been performed by relying on the assumption of a globally 
stationary wind with a smooth density/velocity stratification. 
Consequently, the underlying models are incapable in principle 
of describing the aforementioned features, and the derived 
results (including the verification of the theory) may depend 
strongly on this assumption. 

Theoretical effort to understand the nature and origin 
of these observational findings have generally focused on 
the line-driving mechanism itself; the first linear-stability 
analyses showed the lin e force to be inherently unstable 
JOwocki & Rvbickill 19841 and references therein). Subsequent 
numerical simulations of the non-linear evolution of the line- 
driven flow instability (for a review, see ■OwockL.1994.> . with 
various degrees of approximati on concerning the stabilizing 
diff'use, scattered radiation field ('Ow ocki & Pulsl ll996"l999'). 
have shown that the outer wind (typically, from 1.3 /?* on) de- 
velops extensive structure, consisting of strong reverse shocks 
separating slower, dense shells from high-speed rarefied re- 
gions. Only a very small fraction of material is accelerated to 
high speed and then shocked; for most of the flow the major 
effect is a compression into narrow, dense "clumps", separated 
by large regions of much lower density. 

At first glance, these models appear to be in strong con- 
trast with our assumptions for the "standard model" for wind 
diagnostics based on stationarity and homogeneity, especially 
when viewed with respect to the spatial variation of velocity 
and density. However, when viewed with respect to the mass 
distribution of these quantities, the models are not so very dif- 
ferent (e.g., IOwockietal]ll988t IPuls et alJ[l993al) . Given flie 
intrinsic mass-weighting of spectral formation, and the exten- 
sive temporal and spatial averaging involved, the observational 
properties of such structured models are quite similar to what 
is derived from the "conventional" diagnostics, in an average 
sense. Structured winds also explain 'steady-state' character- 



istics like X-rays, and the black absorption t rough s observed 
in saturat ed UV resonance lines ( iLucvll 19821 Il983l later con- 
firmed bv lPuls et al.ll993al on the basis of hydrodynamical sim- 
ulations). 

Recent time-dependent simula tions have aimed at investi- 
gating two specific problems. First. lRunacres & Owockil 1 120021 
■2005.) have introduced new methods to numerically resolve 
even the outermost wind. In particular, they provide theoretical 
predictions for the radial stratification of the so-called clump- 
ing factor. 



/cl 



< p^ > 
<p >2 



> 1, 



(1) 



where angle brackets denote (temporal) average quantities. For 
self-excited instabilities (e.g., without any photospheric distur- 
bances such as pulsations or sound-waves), they find that, be- 
ginning with an unclumped wind in the lowermost part (/d = 
1), the clumping becomes significant (/^ ^ 4) at wind speeds of 
a few hundreds of km s"', reaches a maximum (/d - 15. . .20), 
and thereafter decays, settling at a factor o f roughly four again. 

On the oflier hand. lPessart & Owockil lEo 03. 2005), build- 
ing on a pilot investigation by Owocki ( 1999), have taken the 
first steps towards including 2-D effects of the radiation field 
into a higher-dimensional hydrodynamical description, to ob- 
tain constraints on the lateral extent of clumping. 

Taken together, and with respect to NLTE modeling and 
spectral analysis, the above scenario has the following ma- 
jor implications, related to radiation field and density/velocity 
structure :- 

1 . X-ray emission arising from the formation and interaction 
of clumps and shocks, in conce rt with an enhanced EUV 
flux (e.g.. iFeldmeier et al.lll997l) . can have a strong influ- 
ence on the ionization/excitation balance in the wind (e.g., 
Pauldrachet al. 2001). 

2. Clumping introduces depth-dependent deviations from a 
smooth density structure, which particularly affects com- 
mon observational mass-loss indicators, such as Hq, emis- 
sion and the IR/radio excess, since these diagnostics di- 
rectly depend on < p^ > (being larger than < p >^). 
Furthermore, the ionization balance becomes modified, pri- 
marily because of the additional < p ^ >-dependence of ra- 
diative recombination rates (see also lBouret et al.l2005l) . 

3. Not only the modified density stratification, but also the 
highly perturbed velocity field can affect the spectral line 
formation, because of its multiple non-monotonic nature. 
This gives rise to modified escape probabilities and multi- 
ple resonance zones for certain frequencies. A major exam- 
ple of such an influence is the formation of black absorp- 
tion troughs in saturated UV-resonance lines (see above). 
Optical lines (e.g. , H„) can also be afl'ected, though to a 
lesser extent (e.g. JPuls et al.ll993bl) . 



Although the potential effects of clumpin 
first discussed some time ago (e.g.. lAbbott et 



ing 
aU 



were 



iLamers & WatersI Il984bt IPuls et alJ Il993bl) . andTave been 
accounted for in the diagnost ics of Wolf-Rayet star s sinc e 
pioneering work by iHillieri (119911) and Schmutzi lll995l) . 
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this problem has been reconsidered by the "OB-star com- 
munity" only recently, mostly because of improvements 
in the diagnostic tools, and particularly the inclusion of 
line-blocking/blanketing in NLTE atmospheric models. 

[Repolust et al. (2004) presented results of a re-analysis 
of Galactic O-stars , previously mo deled using unblanketed 
model atmospheres fP uls et allll996l) . with new line-blanketed 
calculations. As a result of the line blanketing, the derived 
effective temperatures were significantly lower than previ- 
ously found, whereas the modified wind-momentum rates re- 
mained roughly at their former values. Based on this investi- 
gation, and deriving new s pectral-type- rpff an d spectral-type- 
log jg c alibrations (see also Martins et alJ2005l) . lMarkova et alJ 
ll2004 "Paper I") extended this sample considerably and ob- 
tained a "new" empirical wind-mornentum-luminosity rela - 
tionship 1 ( WLR; 'Kudri tzkiet ai]ll995[ iKudritzki & Pui3 l2C 
for Galactic O stars, based on Hq. mass-loss rates. 

A comparison of the "observ ed" wind-moment um rates 
with theoretical predictions from IVink et al.l ( 1200 0!) and in- 
dependent calculations performed by our group (PulsetalJ 
120031) ^. revealed that objects with H[, in emission and those 
with Hq. in absorption form two distinct WLRs. The latter is in 
agreement with theory, whilst the former appears to be located 
in parallel, but above the theoretical relation. This difference 
was interpreted as being a consequence of wind clumping, with 
the contribution of wind emission to the total profile being sig- 
nificantly different for objects with H„ in absorption compared 
to those with H„ in emission (since for the former group only 
contributions from the lowermost wind can be seen, whereas 
for the latter the emission is due to a significantly more ex- 
tended volume). Thus, there is the possibility that for these ob- 
jects one sees directly the effects of a clumped wind, which 
would mimic a higher mass-loss rate (as is most probably the 
case for Wolf-Rayet winds). With this interpretation, the pres- 
ence of clumping in the winds of objects with Hq, in absorption 
is not excluded; owing to the low optical depth, however, one 
simply cannot see it, and coiTesponding mass-loss rates would 
remain unaffected. 

The "actual" mass-loss rates for objects with Hq in emis- 
sion can then be estimated by shifting the observed wind- 
momentum rates onto the theoretical predictions, with a typical 
reduction in M by factors of 2-2.5, corresponding to clumping 
factors of the order of 4-6. 

Though factors between two and three seem reasonable 
when compa red to results from Wol f-Rayet stars (also factors 
of ~3; e.g., Moffat & Robert 1994), there is increasing evi- 
dence that the situation might be even more extreme. From 
an analysis of the ultraviolet Pv resonance doublet (which is 
unsaturated, and can therefore be used as a mass-loss indica- 
tor), MassaeLaL<20Q3) and Fullerton et al. (2004. 2006) con- 
clude that typical O-star mass-loss rates derived from Hq or 



radio emission might overestimate the actual values by factors 
of up to 100 (with a median of 20, if Pv were the dominant io n 
for spectral types be tween 04 to 07; see lFullerton et alJ2006l) . 
iBouret et al.l (120051) . from a combined UV and optical analy- 
sis, obtained factors between 3 and 7, though from only two 
stars. In addition, the latter work suggests that the medium is 
clumped from the wind base on, in strong contrast with typi- 
cal hydrodynamical simulations (see above). If this were true, 
presently accepted mass-loss rates for non-supergiant stars also 
need to be revised, and even the analysis of quasi-photospheric 
lines (i.e., stellar parameters and abundances) might be af- 
fected, since the cores of important lines are formed in the tran- 
sonic region. 

In this paper, we attempt to undertake a first step towards 
a clarification of the present puzzling situation. From a si- 
multaneous analysis of Hq, IR and radio observations, we ob- 
tain constraints on the radial stratification of the clumping fac- 
tor , and test ho w far the results meet the predictions given 
bv lRunacres & Owocki(.200Z.2005.) . Since all these diagnos- 



' The presence of such a relationship is explained by the radiation- 
driven wind theory, namely that the modified wind-momentum rate, 
Mv„(Rt/Roy'^ , should depend almost exclusively on the stellar lumi- 
nosity, LJLq, to some power. 

^ which proved to be almost identical, though the two approaches 
are rather different; see also Kudritzki (2002i) . 



tics depend on < p >, however, we are able to derive only 
relative, not absolute, values, as detailed in Sect. |3 Let us 
point out here that our analysis is based upon the assumption 
of small-scale inhomogeneities redistributing the matter into 
overdense clumps and an (almost) void inter-clump medium, 
in accordance with (but not necessarily related to) the basic 
effects of the line-driven instability. Indeed, the question of 
whether the wind material is predominantly redistributed on 
such small scales and not on larger spati al scales (e.g. , in the 
form of co-rotating intera ction regions; lMullanlll984 Il986t 
ICranmer & Ow ocki 1996) has not yet been resolved, but un- 
explained residuals from the results of our analysis might help 
to clarify this issue. 

Investigations such as we perform here are not new. 
Indeed, a number of si milar studies have b e en presented 
during r ecent years, e.g .. iLeitherer e| al| (Il982l): lAbbott et alJ 
( ^ 1 ^ 8 41 ); JLamers & Leitherer ( 1993); Runacres & Blomme 
Jl996l): feloinine & Runacres. C1997.) : .Scuderi et aL C1998.) : 
iBlomme et al.l (120021 l2003l) . The improvements underpinning 
our study, which hopefully will allow us to obtain more 
conclusive results, are related to the following facts. First, 
and in contrast to earlier work, the uncertainties concerning 
the adopted stellar parameters have been greatly reduced, 
since they have been derived by means of state-of-the-art, 
line-blanketed models. Secondly, we do not derive (different) 
mass-loss rates from the different wavelength domains based 
on a homogeneous wind model, but aim at a unique solution by 
explicitly allowing for clumping as a function of radius, at least 
in a simplified way. Thirdly, we use recent radio observations 
obtained with the Very Large Array (VLA), which, because of 
its gain in sensitivity due to improved performance (mostly 
at 6cm, where the system temperature improved from 60 to 
45 K) allows us to measure the radio fluxes for stars with only 
moderate wind densities, which produce Hq, in absorption. In 
this way we are able to test the above hypothesis concerning 
the differences of Hq mass-loss rates from stars with Hq 
emission and absorption. Lastly, our IR analysis does not 
depend on assumptions used in previ ous standard methods 
exploiting the IR excess (e.g.- lLamers & Waters 1984a) . since 
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Table 1. Sample stars and stellar/reddening parameters as used in this study. Note that radii, mass-loss rates (assuming an un- 
clumped medium) absolute visual magnitudes, Mv, and reddening parameters have been modified with respect to the original val- 
ues (from "ref 1" and "ref2") by a combined V/J/H/K-hand de-reddening procedure (see Sect. 12. 6L using the indicated distances. 
Gravitational accelerations, logg, are "effective" values, i.e., without centrifugal correction, derived from Hy or calibrations; Voo 
is in km s"'; M is in 10"^ Mg/yr; and distances are in kpc. "pt" denotes the Hq, profile type (emission/absorption/intermediate), 
"refl" and "ref2" refer to the sources of the original stellar and magnitude/color/reddening parameters, respectively, where the 
extensions given for reference "1" denote the "preferred" model chosen in Paper I (see text). For f Pup (HD 66811), we provide 
two entries, based on different distances (see Paper I). 



Star 


Sp.Type 


Teff 


logs 


R, 


KHe 


Voo 


pt 


M(opt) 


y8(opt) 


Mv 


E(B-V) 


Rv 


dist 


refl 


ref2 


CygOB2#7 


03If* 


45800 


3.93 


15.0 


0.21 


3080 


e 


10.61 


0.77 


-5.98 


1.77 


3.00 


1.71 


2 


5 


HD 190429 A 


04If+ 


39200 


3.65 


22.7 


0.14 


2400 


e 


16.19 


0.95 


-6.63 


0.47 


3.10 


2.29 


1 


1-0 


HD 15570 


04If+ 


38000 


3.50 


24.0 


0.18 


2600 


e 


17.32 


1.05 


-6.69 


1.00 


3.10 


2.19 


4 


6 


HD66811 


04I(n)f 


39000 


3.60 


29.7 


0.20 


2250 


e 


16.67 


0.90 


-7.23 


0.04 


3.10 


0.73 


3 


1-4 










18.6 








8.26 




-6.23 


0.04 


3.10 


0.46 


3 


1-0 


HD 14947 


05If+ 


37500 


3.45 


26.6 


0.20 


2350 


e 


16.97 


0.95 


-6.90 


0.71 


3.10 


3.52 


3 


1-2 


CygOB2#ll 


05If+ 


36500 


3.62 


23.6 


0.10 


2300 


e 


8.12 


1.03 


-6.67 


1.76 


3.15 


1.71 


2 


5 


CygOB2#8C 


05If 


41800 


3.73 


15.6 


0.13 


2650 


a 


4.28 


0.85 


-5.94 


1.62 


3.00 


1.71 


2 


5 


CygOB2#8A 


05.51(f) 


38200 


3.56 


27.0 


0.14 


2650 


i 


11.26 


0.74 


-6.99 


1.63 


3.00 


1.71 


2 


5 


HD210839 


06I(n)f 


36000 


3.55 


23.3 


0.10 


2250 


e 


7.95 


1.00 


-6.61 


0.49 


3.10 


1.08 


3 


1-2 


HD192639 


07Ib(f) 


35000 


3.45 


18.5 


0.20 


2150 


e 


6.22 


0.90 


-6.07 


0.61 


3.10 


1.82 


3 


1-0 


HD34656 


0711(f) 


34700 


3.50 


25.5 


0.12 


2150 


a 


2.61 


1.09 


-6.79 


0.31 


3.40 


3.20 


1 


1-6 


HD24912 


07.5III(n)((f)) 35000 


3.50 


24.2 


0.15 


2450 


a 


2.45 


0.80 


-6.70 


0.33 


3.10 


0.85 


3 


1-2 


HD203064 


07.5III 


34500 


3.50 


12.4 


0.10 


2550 


a 


0.98 


0.80 


-5.23 


0.23 


3.10 


0.79 


3 


6 


HD36861 


08III((f)) 


33600 


3.56 


14.4 


0.10 


2400 


a 


0.74 


0.80 


-5.52 


0.08 


5.00 


0.50 


1 


1-1 


HD207198 


09Ib/II 


36000 


3.50 


11.6 


0.15 


2150 


a 


1.05 


0.80 


-5.15 


0.58 


2.56 


0.83 


3 


1-1 


HD37043 


09III 


31400 


3.50 


17.9 


0.12 


2300 


a 


1.03 


0.85 


-5.92 


0.04 


5.00 


0.50 


1 


1-1 


HD30614 


09.5Ia 


29000 


3.00 


20.7 


0.10 


1550 


e 


3.07 


1.15 


-6.00 


0.25 


3.10 


0.79 


3 


1-2 


CygOB2#10 


09.51 


29700 


3.23 


30.7 


0.08 


1650 


i 


2.74 


1.05 


-6.95 


1.80 


3.15 


1.71 


2 


5 


HD209975 


09.5Ib 


32000 


3.20 


14.7 


0.10 


2050 


a 


1.11 


0.80 


-5.45 


0.35 


2.76 


0.83 


3 


1-1 


References: l.tMarkova et al. 


(2004), 2. N 


lokiemetal. (2005), 3. Renolust et al., ^2004), 4. .Repolust et ab k2Q.Q.5.). 5., Hanson (i2003.) 
is-Acellaniz et alJJ2004h. 


(distan 


fromlMassev & ThomnsorJ Jl99lhl. 6.lMa 





we calibrate against results from line-blanketed NLTE models. 
(Note that uncertainties in the stellar radii due to distance 
errors cancel out as far as the derived run of the clumping 
factors is concerned, and affects "only" the absolute mass-loss 
rates.) 

The plan of this paper is as follows. In Sect.|2l we describe 
our stellar sample and the observational material used in this 
study. We also comment on some problems related to redden- 
ing. In Sect. 13 we present the methods used to analyze the 
different wavelength regimes, and discuss how we deal with 
a clumped wind medium. Applying these methods, we derive 
constraints for the radial stratification of the clumping factor in 
Sect. |4] and give a discussion and summary of our findings in 
Sects. |5landil 



2. Stellar sample and observational material 

The stellar sample consists of 19 Galactic supergiants/giants, 
covering spectral types 03 to 09.5. These stars have been an- 
alyzed in the optical and, to a large part, (re-)observed by us 
with the VLA. To our knowledge, the only confirmed non- 
thermal radio emitter included in our sample is Cyg OB2#8 A 
(Biegin get all 119891) . which was recently detected as an 
061/05. 5III, colliding-wind binary system by de Becke r et al.l 
t2004) . Somewhat inconsistently, we will use corresponding 



stellar parameters resulting from an analysis assuming a sin- 
gle star. Note also that HD 37043 is listed as an SB 2 binary 
in the recent Galactic O-star catalogue of Mais-Apellaniz et alJ 
ll20n4h . 

Most of the opti cal analyses were perfo rmed by either 
iRepolust et J] ll2004l) or lMokiemet"ai] (12005 ') (Cyg OB2 ob- 
jects), usi ng the NLTE line -blanketed model-atmosphere code 
FASTWiND (Puis et al.U2005l) . For a few stars (those denoted by 
"1" in Tabled column "refl"), stellar parameters have been 
derived from calibrations only, as outlined in Pap er I. At least 
for HP 190429 A, an independent re-analysis bv iBouret et alJ 
(l2005J) . by means of the alternative model-atmosphere code cm- 
FGEN iHillier & Miller, 1998), confirms the corresponding cali- 
bration. Finally, for HD 15570, we use parameters de rived from 
H- & K-hand spectroscopy bv lRepolust et alJ ( l2005l) . 

The parameters adopted in this study are presented in 
Table ^ For those objects which have been analyzed exclu- 
sively in Paper I, and for which more than one choice concern- 
ing distance, reddening or luminosity has been discussed, we 
have used the "preferred" parameter set (Paper I, Table 2), de- 
noted by the corresponding extension in entry "ref 2". Only for 
HD 668 11 (^ Pup) do we provide two entries, referring to its 
"conventional" distance, d - 460 pc (2nd entry), and the as- 
sumption that this star is a runaway star, located at d - 730 pc 
('see Sahu & Blaauw 1993i and Paper I, Sect. 5). Unless stated 
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explicitly, we will use the latter parameter set in our further 
discussion. 

Note that due to minor revisions with respect to redden- 
ing, the stellar radii and H„ m ass-loss rates (rescaled by as- 
suming M/R];^ - const, e.g.. Pu is et al.l 19961) for most objects 
are (slightly) different from the original sources. In Sect 12.61 
we will discuss why these revisions were necessary, and how 
they have been obtained. 

2.1. Variability of the diagnostics used 

Before we discuss the observations obtained in the individ- 
ual bands (H„, IR and radio), let us first give some important 
comments on the variability of the different diagnostics. Stellar 
winds are known to be variable on different timescales and in 
all wavelength ranges in which they are observed. Thus, the use 
of non-simultaneous measurements, as in our analysis, can be 
an issue. 

Regarding Hq,, line profile variations in early-type stars 
have be en observed for years . Since the first extensive sur- 
veys by iRosendhall ( Il973albl) . a large number of investiga- 
tions have been conducted to establish the properties of the H„ 
variability and also its origin (e.g., Ebbets 1982; Sc uderi et al. 
ITggilKaufer et al Jl 996HKaner et al Jl 997HMorel";t al. 2004). 
Although the variations in the Hq, profile in some cases look 
very dramatic, they indicate, when interpreted in terms of a 
variable mass-loss rate, only moderate changes in M, usually 
not exceeding the uncertainties on the corresponding estimates. 
Recen tly, for a sample of 15 0-type supergiants, Markova et al. 
1I2OO5I) constrained the M variability to about +4% of the cor- 
responding mean value for stars with stronger winds, and to 
about +16% for stars with weaker winds. These estimates are 
in r emarkably good agreement with th ose from previous stud- 
ies (lEbb etsHl982HScuderi et allll992l) who report variations in 
M of about 10 to 30%. 

In the case of IR and radio continua, and assuming the 
emission to be thermal, the timescales of variability (due to 
variations of micro- or macro- structure, i.e., of the local den- 
sity or mass-loss rate^) can differ by orders of magnitude in the 
two wavelength regimes. Considering variations in M, the tran- 
sit time of a front would be of the order of hours in the near-IR 
forming region (given typical sizes of the emitting region and 
velocity of the expanding material), and as much as months, 
or even years, in the radio domain. This implies that whilst the 
IR emission would display short-term variability, following the 
mass-loss rate variations very closely, variations in the radio 
would be averaged out if they occurred on timescales much 
shorter than the transit time. 

Different considerations apply when the variability is of 
non-thermal origin. In this case, only the radio emission is af- 
fected. The proces s responsible i s usually cited as being syn- 
chrotron emission (White, 1985), most probably produced in 
coUiding-wind binaries ( Van Loo et al., 2006) . The main ob- 
served characteristics are variability over timescales of up to 
months, and a power-law spectrum increasing with wavelength 



^ Note that variations in the ionization can also induce temporal 
variability, e.g.. Panagia 1 1991) . 



and with a variable spectral index ( iBieging et all 119891) . In 
such a case, which is met at least by one of our objects, 
Cyg OB2#8 A, the measured radio-flux(es) can still be used as 
an upper limit of the thermal free-free emission, by analyzing 
the lowermost flux measured at the shortest radio wavelength. 

Regarding the amplitude of variability, no clear evidence of 
IR continuum variability has been reported up to now. Amongst 
the IR observations we have obtained from the lit erature, there 
are so me studies (e.g.. ICastor & Simonlll983l or lAbbott et alJ 
Il984l) with data sampled on timescales ranging from a few 
hours up to a few months, but no variation of the observed 
IR fluxes above the errors was reported. If, on the other hand, 
we compare sets of measurements of the same object, ob- 
tained by different authors with different instruments, we do 
observe differences in the measured fluxes, more likely related 
to calibration problems than to genuine IR variability (see also 
Sect. l2ni . 

With respect to radio emission, there are several pieces of 
evidence for variability, both in the observed fluxes and in the 
spectral index. Again, we have to distinguish between thermal 
and non-thermal emission. In th e case of non-therrn al origin, 
variability is always present (e.g., Bieging et alJl9 89). This has 
to be accounted for whenever we have no clear indication about 
the thermal origin of the observed emission, but where we do 
see variations. Of our targets, in addition to #8A, this might 
be a problem only for HD 190429A (and for HD 34656 and 
HD 37043 for oflier reasons). 

Among thermal emitters, on the other hand, the situation 
is less clear. There are very few studies which have observed 
one object several times and at more tha n one frequency. In the 
sample studied by Bieging et al. (1989), two from six definite 
thermal emitters showed variability, both of which are B super- 
giants (CygOB2#12 and ^ Sco). In these cases, the flux vari- 
ation reached values of up to 70%. Interpreted in terms o f M, 
this w ould mean a change of 50% (see Eq.|5J. In Scuderi et alJ 
(Il998l) . one out of six objects (again CygOB2#12) showed 
variability whilst ha ving a spe ctral index compatible with ther- 
mal emission. Blomme et alJ (2002) studied the variability of 
e Ori (BOIa) and found no evidence for variability, both on 
shorter and longer timescales. The best-studied object with re- 
gard to therr nal radio variability is ( Pup (04If-i-), as a result of 
the work bv lBlomme et al.l (l2003l) . who investigated both new 
and various archival data. Again, short-term variability could 
be ruled out, and long-term variability (with observations be- 
ginning in 1978) appeared to be low or even negligible. 

The major hypothesis underlying our present investigation (be- 
ing in agreement with most other investigations performed thus 
far) is that the clumping properties of a specific wind are con- 
trolled by small-scale structures. Further comments on this hy- 
pothesis, in connection with the outcome of this analysis, will 
be given in Sect. |6] If related to any intrinsic wind property 
(e.g., the instability of radiative line-driving, even if externally 
triggered by short wavelength/short period modulations), the 
derived clumping properties should be (almost) independent of 
time, as long as the major wind characteristics remain largely 
constant. Accounting for the observational facts above, this as- 
sumption seems to be reasonable, and justifies our approach 
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Table 2. VLA radio fluxes (in juJy), with 1-cr errors in brackets. Data without superscripts are new observations (see Sect. 12. 3> . 
whilst data with superscripts correspond to either (a) unpublished rneasure ments by Scuderi et al. or liter ature values used to 
compl ement our database: (b) IScuderiet a l. ( 1998); (c) Bieaina et al. ( 1989); (d) 3.6 cm observations from Lamers & Leitheren 
lll993h (concerning HD 15570, see text); (e) .Blomme et aL (.2003. including 20 cm data for ^ Pup, at 760+90 //Jy). 
Also indicated are the adopted IR to mm fluxes and the sources from which they have been drawn (see foot of table). Data 
denoted by "own" refer to JHKLM-hand observations performed by OGT at the Crimean 1 .25 m telesco pe (see Sec t. 12. 4> . scuba 



data (at 1.35 mm) were obtained by AWE and IDH (see Sect.|2.5|l 


. and SCUBA 0.85 mm data are from,Blomme et ali t2003.). 


Star 


4.86 GHz 


8.46 GHz 


14.94 GHz 


43.34 GHz 


IR- and mm- 


references 




(6 cm) 


(3.5 cm) 


(2 cm) 


(0.7 cm) 


bands used 


(IR and mm) 


CygOB2#7 


<112 


<100 






HKLMN 


1,14 


CygOB2#8A 


<540" 

1000(200)'- 

800(100)'' 

700(100)'- 

400(100)'- 


920(70)" 


500(200)'- 




JHKLMNQ 


1,5,14,19,20 


CygOB2#8C 


<200'- 








HKLMN 


14 


CygOB2#10 


134(29) 


155(26) 


300(100) 




JHKLMN 


5,14,19 


CygOB2#ll 


182(33) 


228(28) 


<400 




JHKLMN 


5,14 


HD 14947 


<110 
<90" 


<135 
90(30)" 
120(30)" 
110(30)* 


<700 




JHKLMN 


2,5,15,own 


HD 15570 


100(40)" 


220(40)" 

125(25)'' 






JHKLMNQ,l.35mm 


1,5,8,1 l,15,18,scuBA 


HD24912 


<200 


<120 


<390 


<840 


JHKLMNJRAS 


3,5,7,16 


HD 30614 


230(50)* 


440(40)* 


650(100)* 




JHKLMN 


5,7,own 


HD 34656 


<132 


119(24) 


<510 




JHKL 


17, own 


HD 36861 


<112 


<90 


<1000 




JHKLMN 


2,5,7 


HD 37043 


203(38) 


<90 
46(15)'' 


<330 




JHKLMN 


4,5,16,21,22 


HD66811 


1640(70)' 


2380(90)'' 


2900(300)" 




i//A:LM,IRAS,0.85mm,1.3mm 6,9,10,12,13,22,23,24 




1490(110)" 












HD 190429 A 


250(37) 


199(36) 
280(30)* 


<420 


<540 


JHKLM 


5,20,own 


HD 192639 




<90" 






JHKLM 


5, 15, own 


HD 203064 


114(27) 


126(20) 


<330 




JHKLMJRAS 


3, 5, own 


HD 207 198 


105(25) 


101(21) 


249(82) 




JHKLMJRAS 


3, own 


HD 209975 


165(36) 


184(28) 


422(120) 




JHKLMJRAS 


3, own 


HD 210839 


238(34) 


428(26) 


465(120) 


790(190) 


J HKLMNQJRAS, 1. 35 mm l,2,3,5,14,15,own,scuBA 


References for IR and mm data: 1. Abbott et al. ( 1984), 


2. Barlow & ( 

Gehrz et al. I 

1.1 1966b), 13. 

Sagar & Yu ( 1 


:;ohen (1977), 3. Beichman et all ll988h. 4. Breser et al. (1981). 5. 


Castor & SimoiJ h983h. 6.lDachs & Wams 


teker ( 1982), 7. 
12.Johnson et a 


1974). 8. Guetter & Vrbt 


h989ft. 9.ljohnson & Borsman ( 1963). 10. 


.lohnsorlJl964h. 


ll.|johnsonetalJh966a1. 


Lamers et al. ( 198^), 14.lLeitherer et aLtl98^). 15.[Leitherer & Wold 


il984h. 16.lNev 


etalJJl973h. 17 
mBreda|jl980^ 


.|PolcaroetalJil990). 18. 


990). 19. Sneden et alJJl978h. 20.lTaDiijll98ll). 21.lThe et alJJl98d). 


22.|Whittet&v. 


23.|Leitherer& RobertH1991D, 24.|Blomme et alJ^2003^. 





of using observational diagnostics from different epochs. Note 
also that the observed X-ray variability (where the X-rays are 
thought to arise mostly from clump-clump collisions) is low 
as well (Berghoferet al.. 1996) . due to the cancellation effects 
of the large number of participating clumps being accelerated 
in laterally independent cones (of not too large angular extent, 
[Peldmeieretal. 1997). 



If we had analyzed only one object, the derived results 
might be considered as spurious, of course. However, due to the 
significant size of our sample, any global property (if present) 
should become visible. Let us already mention here that our 
findings, on average, indicate rather similar behaviour for sim- 
ilar objects, and thus we are confident that these results remain 



largely unaff'ected by issues related to strong temporal variabil- 
ity. 



2.2. Ha observations 



For our analysis of H„ by means of clumped wind models, 
we have used the same observational material as described 
in Paper I, i.e., H„ spectra obtained at the Coude spectro- 
graph of the 2m RCC telescope at the National Astronomical 
Observatory, Bulgaria, with a typical resolution of 15 000. For 
further information concerning technical details and reduction, 
see Paper I and references therein. 
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2.3. Radio observations 

New radio observations for 13 stars have been carried out at the 
VLA (in CnB and C configuration), in several sessions between 
February and April 2004, for a total of about 36 hours. Some 
of these stars were already known to be radio emitters, but for 
many of them only upper limits for their radio emission were 
available. Exploiting the gain in sensitivity of the VLA, and 
guided by the requirement of using consistent data at all radio 
frequencies for our analysis, we decided to (re-)observe them. 
Note particularly that it was possible to observe not only stars 
with strong winds, but also those with weaker winds (i.e., with 
Hq. in absorption). 

The journal of observations is given in Table lA.ll with 
dates, observing frequencies, time on targets, calibrators for 
flux-density bootstrapping and VLA configuration. The obser- 
vations were performed with a total bandwidth of 100 MHz at 
all frequencies. The target stars were observed for several scans 
of about 10 minutes, interleaved with a proper phase calibrator 
at 4.86, 8.46, and 14.94 GHz. A faster switching between the 
target star and the phase calibrator was used at 43 GHz, to re- 
move the rapid phase fluctuations introduced between the an- 
tenna elements by the troposphere at this frequency. The data 
at 15 and 43 GHz have been corrected for atmospheric opac- 
ity using a combination of a seasonal model and the surface 
weather conditions during the experiment. The Astronomical 
Image Processing System (AIPS) developed by the National 
Radio Astronomy Observatory (NRAO) was used for editing, 
calibrating and imaging the data. 

Table |2l (left part) displays the corresponding fluxes, to- 
gether with data from other sources (in most cases, at 2, 3.5 
and 6 cm) used to complement our sample. For four objects, 
partly overlapping with the 13 stars mentioned above, we re- 
lied on data (denoted by superscript "a" in Table |2} derived by 
Scuderi et al. (in preparation), using the VLA in B, BnC and 
C configuration, in different sessions between 1998 and 1999. 
The reduction and analysis of these data has been performed in 
a similar way as outlined above for our new observations. The 
quoted flux limits for both datasets refer to 3 times the RMS 
noise on the images, whereas the errors are 1-cr eiTors, again 
measured on the images. Note that at these low flux densities 
the contribution of errors introduced by the calibration proce- 
dure is negligible. 

For the rem aining stars, li t eratur e val ues have been used, in 
particular from lScuderi et alJ (Il998l) and'Bieaing et al.' ('1989^, 
togeth er with 3.6 cm observations from Lamers & Leithereg 
(I1993I) : these are denoted by superscripts "b", "c" and "d", re- 
spectively. Note that the indicated 3.6 cm flux for HD 15570 
deviates from the "o r iginal" value of 1 10 + 30 fiiy provided by 
iLamers & Leithereii lll993ll . as a result of a recent recahbra- 
tion of the original Howarth & Brown VLA data, performed 
by IDH. For i f Pup , finally, we used the data obtained by 
iBlomme et al.l (l2003[ denoted by superscript "e"), at 3.6 and 
6 cm (Australia Telescope Compact Array, ATCA) and 20 cm 
(VLA), in combination with the 2 cm data from Bieging et al. 
(recalibrated, see Blomme et al.) 

For those objects which have been observed both by us 
and by others, or where multiple observations have been ob- 



Table 3. Near IR magnitudes and eiTors (last two digits) for 
program stars as observed with the Crimean 1.25 m telescope. 



Star 


JD 

(2453+) 


J 


H 


K 


L 

M 


HD 14947 


067.204 


7.17 02 


6.95 02 


6.88 01 


6.85 04 




307.456 


7.10 01 


6.98 01 


6.85 04 


6.67 08 


HD 30614 


073.238 


4.32 02 


4.25 02 


4.20 02 


4.20 01 




100.224 


4.30 01 


4.26 01 


4.27 01 


4.23 01 


HD 34656 


072.356 


6.67 01 


6.71 01 


6.6101 


6.60 04 




100.244 


6.69 02 


6.71 01 


6.69 01 


6.65 03 


HD 190429 A 


216.414 


6.28 01 


6.12 01 


6.1401 


6.13 04 




225.439 


6.18 01 


6.01 01 


6.1901 


5.98 08 


HD 192639 


216.439 


6.45 01 


6.24 01 


6.22 01 


6.26 04 




307.254 


6.44 01 


6.23 01 


6.17 01 


6.24 04 


HD 203064 


223.459 


5.17 01 


5.12 01 


5.13 01 


5.13 03 
4.98 10 




311.301 


5.19 01 


5.17 01 


5.17 00 


5.02 02 
5.02 05 


HD 207198 


223.496 


5.5101 


5.35 02 


5.39 01 


5.37 04 
5.58 10 




309.167 


5.48 02 


5.42 01 


5.45 01 


5.58 03 
5.57 20 


HD 209975 


223.535 


5.01 01 


4.97 01 


5.00 01 


5.12 05 
5.00 07 


HD 210839 


223.567 


4.62 01 


4.52 01 


4.54 01 


4.57 02 
4.44 05 




309.193 


4.61 01 


4.51 01 


4.58 01 


4.62 02 
4.37 06 



tained (particularly for the non-thermal emitter Cyg OB2#8A), 
we have added these values to our database. In almost all cases, 
the different values are consistent with each other, especially 
for the weaker radio sources w hen co mparing with the upper 
limits derived bv .Bieg in g etalJ lll989l) . 



2.4. IR observations 

In the right part of Table |2l we have summarized the IR data 
used, which are to a large part drawn from the literature. For 
a few objects, IRAS data for 12 , 25, 60 and 100 /im are also 
available CBeichman et all Il988l) . unfortunately mostly as up- 
per limits for A >25 /im. For ^ Pup (HD 6681 1), however, ac- 
tual values_^em'esen^t all but the last wavelength (100 fim); 
see lLamers et alJ ( Il984 . 

For nine objects (denoted by "own" in the "references" col- 
umn of Table 13, new JHKLM fluxes (see Table |3} have been 
obtained at the 1.25 m telescope of the Crimean Station of 
the Sternberg Astronomical Institute (Cassegrain focus, with 
an exit aperture of 12"), using a photometer with an InSb de- 
tector cooled with liquid nitrogen. Appropriate stars from the 
Johnson catalog (Johnson et al.. 19 66b) were selected and used 
as photometric standards. Where necessary, the HLM magni- 
tudes of the standards have been estimated from their spectral 
types using relations from iKoorneef ( 1983) . 
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2.4.1 . Absolute flux calibration 

In order to convert the various IR magnitudes from the litera- 
ture and our own observations into meaningful (i.e., internally 
consistent) physical units, we have to perform an adequate ab- 
solute flux calibration. For such a purpose, at least three differ- 
ent methods can be applied: - 

1 . Calibration by means of the solar absolute flux, using anal- 
ogous stars. 

2. Direct comparison of the observed Vega flux with a black- 
body. 

3. Extrapolation of the visual absolute flux calibration of 
Vega, using suitable model atmospheres. 

Although the first two methods are more precise, the latter one 
provides the opportunity to interpolate in wavelength, allowing 
the derivation of different sets of IR-band Vega fluxes for var- 
ious photometric systems. Thus, such an approach is advanta- 
geous in the case encountered here (observational datasets ob- 
tained in different photometric systems), and we have elected 
to follow this strategy. 

Atmospheric model for Vega. To this end, we used the lat- 
est Kurucz models^ to derive a set of absolute IR fluxes for 
Vega in a given photometric system, by convolving the model 
flux distribution (normalized to the Vega absolute flux at a 
specific wavelength; see below) with the corresponding filter 
transmission functions. In particular, we used a model with 
Jeff = 9550 K, log^ = 3.95, [M/H] = -0.5 and Vt = 2.0 
km s~' llCaste lli & Kurucz, 1994). In order to account for the 
possibility that the metallicity of Vega might differ from that 
adopted by us, an alt ernative model with [M/H] = -1.0 (cf. 
lOarcia-Gil et alJ2005l) was used to check for the influence of a 
different metallicity on the derived calibration. At least for the 
Johnson photometric system, the differences in the correspond- 
ing fluxes turned out to he always below 1%. 

Visual flux calibration. The most commonly used visu al flux 
calibra tion for Vega is based on the com pilation byjHavesI 
111 9851) . which has since been questioned bv lMegessieii(ll995l) . 
who recommends a value being 0.6% larger than the value pro- 
vided by Hayes (3540 Jy), and equals 3560 Jy (i.e., 3.46 • 10"'' 
erg cm"^ s"' A"') at /I = 5556 A. This value has been used 
when normalizing the Kurucz model fluxes to the monochro- 
matic flux at /I = 5556 A. Since the standard error of the 
Megessier calibration is about one percent, this error is also 
inherent in our absolute flux distribution. 

Vega V-band magnitude. Th e available V-band magnitu des 
of Vega range from 0?026 JBohlin & GiUilandi l2004h to 
0™035 (CoHna & BohUn[ fl994l) . while in the present investiga- 
tion we adopt V = 0™03 mag in agreement with Johnson et al. 
(Il966bl) . With this value, the monochromatic flux for a Vega- 
like star at the effective wavelength of the V filter is F55oo(mv 
= O-^O) = 3693 Jy. 



Filter transmission functions. To calculate the absolute fluxes 
of Vega in a given photometric system, we have to know the 
corresponding filter transmission functions, for each band of 
this system. In those cases where such functions were explicitly 
available we used them, while for the rest (including our own 
IR data) we used trapezoidal transmission curves based on the 
published effective wavelength and FWHM of the filters.^ The 
use of trapezoidal instead of actual response functions might, 
of course, lead to some error in the derived absolute fluxes. 
Indeed, in the particular case of the ESO filt er system, this 
error was estimated to be less than 5% f Schwarz & Melnickt 
I1993I) . with typical values of about 2% systematically larger 
fluxes from the trapezoidal approximation. 

Vega IR magnitudes. To convert stellar magnitudes into ab- 
solute fluxes using Vega as a standard, the magnitudes of Vega 
in the different filters for the various photometric system have 
to be known. In our case, these data have been taken from the 
corresponding literature, and the errors inherent to these mea- 
surements are usually very small. 

Finally, let us mention that we are aware of the problem 
that the use of (simplified) model atmospheres for calculating 
the IR flux distribution of Vega might le ad to s ome uncertain- 
ties, as discussed bv iBoh lin & Gilliland] ( 120041) (e.g., the pos- 
sibility that Vega is a pole-on rapid rotator, Gulliver et al 1994; 
Peter son et al...2004) . Note, however, that .Tokunaga & Vaccai 
(20051) have recently shown that the near-IR (1 to 5 //m) ab- 
solute flux d ensities of V ega derived by means of atmospheric 
models (e.g., ICohen et al. 1992) and by means of direct mea- 
surements (e. g.. lMegessieii 1 9951) are actually indistinguishable 
within the corresponding uncertainties, which, in these specific 
cases, are of the order of 1.45% and 2%, respectively. 

On the other hand, given the fact that Vega has a dust and 

r ir"— — I 

gas disk (IWilner et a l.. 2002) which produces an IR excess, one 
cannot exclude the possibility that a flux calibration based on 
a comparison of Vega observed magnitudes and model fluxes 
might lead to systematic errors, at least for A > 5 ;um, as dis- 
cussed also by .Megessier. C1995.) . There are 12 stars in our 
sample for which we have ground based mid-IR photometry 
obtained in the A^- and Q-bands. In the case that Vega indeed 
displays a mid-IR flux excess (as compared to the models), one 
might expect that the observed fluxes of our targets (based on 
this calibration) are somewhat underestimated in these bands. 
Such a systematic error can be easily detected, however, and we 
shall keep this possibility in mind when performing our analy- 
sis. 



2.5. Mm observations 

For three objects, we were also able to use 1.3/1.35 mm 
fluxes, acquired either with the Swedish ESO Submillimeter 
Telescope (SEST) at La Silla (^ Pup; see iLeitherer & Robert! 
■ 1991.) or wit h the Submillim etre Common User Bolometer 
Array (scuba: iHolland et alJl9 99) at the James Clerk Maxwell 



^ For more detailed information about the shape of the fil- 
ter transmission functions used to convert the literature data, see 



from.http : //kurucz . harvard . edu/stars/vega| 



iRunacres & Blomme 1 1996!. their Table 3). 
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Table 4. 1.35 mm fluxes and errors for program stars observed 

with SCUBA. 



Star 


date of obs. 


integration (s) 


flux (mJy) 


CygOB2#8A 


May 7, 1998 


3600 


-2.50 ±5.95 


HD 15570 


Jul 3, 1998 


2160 


4.76 + 2.43 


HD2I0839 


May 4, 1998 


4500 


4.25 ± 1.92 




Jun 1, 1998 


2340 


8.87 + 3.58 



Telescope (HD 15570 and HD 210839). (For CygOB2#8A, 
which was also observed with scuba, only badly defined upper 
limits were obtained.) 

The SCUBA observations were obtained in the instrument's 
photometry mode (the standard mode employed for point-like 
sources), using the single, 1.35 mm photometric pixel, located 
at the outer edge of the long-wavelength (LW) array. The data 
were acquired in service mode over the period May-July, 1998. 
Table l^lists the observation dates, integration times and mea- 
sured fluxes. Data reduction was performed using the scuba 
User Reduction Facility (surf; Jenness & Lightfoot 2000). 

Additi onal 0.85 mm scuba data have been taken from the 
literature ( iBlomme et alll2003l) . again for f Pup. 

2.6. De-reddening and stellar radii 

Since we are aiming for a combined optical/IR/radio study, all 
parameters used have to be consistent in order to allow for a 
meaningful analysis of the observed fluxes, in particular the 
excesses caused by the (clumped) wind alone. To compare the 
observed wifli the theoretical fluxes, we have (/) to de-redden 
the observed fluxes and (//) to derive a consistent stellar radius 
for a given distance d (or vice versa, see below), which has 
been drawn from the literature cited or recalculated from the 
assumed value of My (for models "1-2" in column "ref2" of 
TableU). 

For this purpose, we have used our (simplified) model as 
described in Sect. l3.2l to synthesize theoretical VJHK fluxes. *■ 
Note that this model has been calibrated to reproduce the cor- 
responding predictions obtained from a large OB-model grid 
calculated by fastwind (Pui s et all 120051) . 

By comparing the observed IR fluxes (from the various 
sources given in Table |3 with the theoretical predictions, we 
derive "empirical" values for the color excess E(B-V) and/or 
the extinction ratio Ry, by requiring the ratio between de- 
reddened observed (plus/minus error) and distance-diluted the- 
oretical fluxes to be constant within the V- to ^-bands. F or this 
purpos e, we adopt the reddening law provided bv Cardelli et al.l 
(Il989l) . Visual fluxes have been calculated using V -magnitudes 
from Paper I or from Mais-Apel laniz et alJ ll2004t . 

In a second step, we adapt the stellar radius (for a given 
distance) in such a way that the mean ratio becomes unity. This 
procedure ensures the correct ratio between radius and dis- 
tance, i.e., angular diameter, which is the only quantity which 
can be specified from a comparison between synthetic and ob- 
served fluxes. Of course, we could have also chosen to modify 
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^ Only near-IR fluxes were used to ensure that the flux excess due 
to the wind remains low, i.e., rather unaffected by clumping. 



Fig. 1. De-reddening procedure, for the example of 
CygOB2#8A. Displayed is the ratio of distance-diluted, 
theoretical fluxes and de-reddened, observed VJHK fluxes, 
as a function of wavelength, with bars accounting for the ob- 
servational errors. Grey entries correspond to "fit-parameters" 
E(B-V) = 1.9, Rv = 3.0 and R, = 24 Rq. Obviously, the extinc- 
tion is too large: the ratio of theoretical to de-reddened fluxes is 
much smaller at shorter than at larger wavelengths (extinction 
decreasing with wavelength). Moreover, the assumed radius 
is too small, since the mean flux ratio (indicated by a dotted 
line) is well below unity (too small an angular diameter). The 
black entries show our final solution, for E(B-V) - 1.63, Rv- 
3.0 and Rt- 27 Rq. Any curvature has vanished, the optical 
flux corresponds to the mean, and the mean ratio itself (again 
indicated by a dotted line) is located at unity. 



the distance for a given radius; however, in order to be consis- 
tent with previous mass-loss estimates from radio observations, 
which rely on certain distances, we have followed the former 
approach. Fig. [Ogives an impression of this procedure, for the 
example of Cyg OB2#8 A. 

Note that in parallel with re-defining the stellar radius, the 
mass-loss rate used to calculate the V- to /T-band model fluxes 
has to be modified as well, because the latter depend on the as- 
sumed value of M (see below). Since we do not know the actual 
mass-loss rate in advance, we follow a simplified approach and 
use a value equal or related to the Hq. mass-loss rate provided 
by previous investigations (entry "refl"). This mass-loss rate, 
however, had been derived for a certain stellar radius, which 
we claim to improve by our procedure. Consequently, we also 
have to modify our "input value" of M, to maintain the Hq. fit- 
quality of the former investigations. As outlined already above, 
this can be obtained by keeping the ratio Q' = M/R];^ constant. 

This scaling has a further advantage, namely that not only 
Hq, but all p^-dependent diagnostics (i.e., H^^ profile shape, IR 
and radio fluxes) and the finally derived run of the clumping 
factor remain almost unaff'ected if a dififerent radius or distance 
(though identical angular diameter) are chosen. 

This notion follows from the fact that the Hq, profile shape 
depends on Q' alone (for given Voo and assuming that the NLTE 
departure coeflicients do not vary), and that the IR and radio 
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optical depths scale with this quantity as well, whereas the cor- 
responding fluxes are additionally diluted by {Ri,/d)^. As an 
example, remember that under certain conditions (see Sect l3.3> 
the radio fluxes scale according to 






(2) 



In other words, as long as Q' and the angular diame- 
ter ("measured" from aligning synthetic and observed, de- 
reddened fluxes; see Fig. ^ remain conserved, almost all fur- 
ther results become independent of the individual choice of 
R, or d, and a translation of our results to different assump- 
tions, e.g., due to future improvements concerning distance 
measurements (gaia), becomes easily possible. The only quan- 
tities which depend directly on these values are the mass-loss 
and wind-momentum rate (e.g.. Paper I), which are of minor 
importance regarding the objectives of this paper. 

One problem inherent to our approach is the fact that the 
derivation of reddening parameters and /?, requires an a pri- 
ori knowledge of M (and clumping properties), since, as stated 
already above, the model fluxes depend on this quantity. 

First note that the flux excess increases as a function of M. 
Consequently, the average slope of the model fluxes decreases, 
which affects our de-reddening procedure (operating in the V- 
to A'-band). This dependence, however, is only moderate, due 
to the rather low excess in this wavelength region for typical 
OB-star winds. Moreover, it is predicted correctly by our mod- 
els if Q' is of the correct order. 

The absolute flux level in the optical and near IR, on the 
other hand, is much more affected by our choice of M, thus 
influencing our derivation of R,. For identical stellar param- 
eters, the V-flux is a (monotonically) decreasing function of 
MJ To a large extent, this behaviour is induced by a decreas- 
ing source function at bf-continuum formation depth, related 
to the decrease in electron temperature (at t(M) w 2/3) when 
M is increasing, and increasing electron scattering. Both ef- 
fects apply to blanketed and unblanketed models; the "only" 
difference concerns the absolute flux level at optical and (N)IR 
bands, which is larger for blanketed models, due to flux- 
conservation arguments (compensation of the blocked (E)UV 
radiation field). 

Since a precise knowledge of the "real" wind density and 
the near-photospheric clumping properties is not possible at 
this stage, only an iteration cycle exploiting the results of our 
following mass-loss/clumping analysis could solve the prob- 
lem "exactly". 

In order to avoid such a cycle, we follow a simplified ap- 
proach, in accordance with our findings from Paper I and antic- 
ipating our results from Sect.0](cf. column "ratio" in Table^}. 
To calculate the theoretical fluxes required for our de-reddening 
procedure, for objects with H„ in absorption we have used the 
actual, 2'-scaled, H„ mass-loss rate, whereas for objects with 



' More precisely: for those wavelength bands where the wind is not 
optically thick, i.e., where the fluxes depend on both the photospheric 
radiation and the wind absorption/emission, there is an additional de- 
pendence on the wind density, oc MIR\, which scales somewhat dif- 
ferently than Q' . 



Hq, in emission we have reduced the corresponding value by a 
factor of 0.48. This approach is based on our hypothesis that 
the lowermost wind is unclumped (see Sect. 13. 4t . and that the 
previously derived Wa mass-loss rates for objects with Hq, in 
emission are contaminated by clumping, with average clump- 
ing factors of the order of ( -^ 1 . 

From the almost perfect agreement of the theoretical V- 
io-K fluxes with the observations for our final, clumped mod- 
els, this assumption seems to be fairly justified. In any case 
(i.e., even if the lowermost wind were to be clumped as well), 
the most important quantity is the effective mass-loss rate (i.e., 
the actual, unclumped M times square root of local clumping 
factor), so any reasonable eiTor regarding this quantity would 
barely affect the corresponding theoretical fluxes and thus our 
de-reddening procedure. 

We will now comment, where appropriate, on the results of 
our procedure for a few individual objects. For the majority of 
stars, only small modifications of the E(B-V) values resulting 
from optical photometry, {B - V), and intrinsic colors, {B - V)o, 
were necessary, while keeping the total-to-selective extinction 
ratio, Ry, at its "normal" value of 3.1, or at a value suggested 
from other investig ations. The iii trinsic colors used here have 
been adapted from Wegneil lll994l) . particularly because of their 
extension towards hotter spectral types. However, since this 
calibration deviates considerably from the widely used alterna- 
tive provided bv .Fitzgerald. 1.1970.) at the cool end (-0.24 mag 
vs. -0.28 mag for 09.5 supergiants), we adopt, as a compro- 
mise, only values < -0.27, and -0.27 if Wegner's calibration 
exceeds this threshold. 

Concerning the CygOB2 stars, for three objects (#7, #8A 
and #8C), our procedure result s in rathe r simila r reddening pa- 
rameters to those presented bv iHansonl ( 12003 1 based on UBV 
photometry by Massey, priv. comm., and IR-photometry from 
2MASS). Only for stars #10 and #11 did we find larger dis- 
crepancies, which were corrected for by using Ry - 3.15 in- 
stead of Ry = 3.0, as suggested by Hanson and previotis wor k 
in the optical (iMassev & Thom pson 1991; Torres et alJll99ll) . 
Note that "our" value is consistent with the values provided by 
Patriarchi et al. (2003, see below): Ry = 3.17 and 3.18, respec- 
tively. 

In disagreement with the work by Hanson, however, we still 
used the canonical distai ice of d - 1.71 kpc for the C yg OB2 
stars, as determined by iMassev & Thompson! (1199 11) . In our 
opinion, the alternative, lower value(s) claimed by Hanson 
would result in too low luminosities. Most probably, though, 
the "real" distance is smaller than the value used here. As 
pointed out, this would imply "only" a down-scaling of radii 
and mass-loss rates, and would not affect our conclusions con- 
cerning the wind clumping. 

The only other objects worthy of closer inspection are those 
with extinction ratios Ry ^ 3.1 (cf. Paper I). Unfortunately, 
the recent catalogue of Ry values for Galactic O-stars by 
IPatriarchi et alJ (2003) covers only few stars in our sample (in 
particular, CygOB2#8A, #10, #11 and HD 34656), such that 
a comparison is not possible for the majority of objects. Due 
to the degeneracy of E(B-V) and Ry (different combinations 
can result in rather similar extinction laws, if the considered 
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Fig. 2. Differences between derived color excess, E(B-V), and 
corresponding literature value, {B - V) - (B - V)o (Table ^ 
entry "ref 2"; for intrinsic colors see text), as a function of 
Teff. Asterisks denote supergiants, and crosses bright giants 
and giants, respectively. The mean deviation for supergiants is 
-0.004 + 0.016 mag, and for I.e. 11/111 stars -0.01 + 0.023 
mag. 



wavelength range is not too large), we applied the following 
philosophy: in those cases with peculiar extinction ratios, Ry 
^3.1, (obtained from Paper 1 and references therein), we firstly 
checked, by using these values, whether the derived color ex- 
cess is consistent (within small errors; see below) with mea- 
sured and intrinsic colors in the optical. If so, we adopted these 
values here. In this way, we confirmed the values of Ry - 5.0, 
5.0 and 2.76 for the stars HD 36861, HD 37043 (both belonging 
to Ori OB 1) and HD 209975, respectively. 

For HD 207 198, with a literature value of Ry =2.76, the de- 
rived color excess would have fallen 0.04 mag below the "op- 
tical" value (a deviation which we considered to be too large if 
Ry ^3.1 anyway). Therefore, for this object, we kept the opti- 
cal E(B-V) value and fitted Ry using our procedure, resulting 
in Ry = 2.56. This star is the only one for which our proce- 
dure showed significant deviations from previous work. Given 
the difficulties in deriving reliable Ry values, however, we con- 
sider this deviation as not too troublesome. 

For the last object in this group, HD 34656, we could check 
for the consistency of our results with the work by Patriarchi et 
al. By keeping Ry - 3.1, as suggested in Paper 1, the derived 
E(B-V) would lie 0.03 mag above the "optical" value, which, 
compared to the other objects (see below), is rather large. On 
the other hand, by keeping our value of E(B-V), we derived Ry 
= 3.4, which is consistent with the value claimed by Partriarchi 
et al. (Ry = 3.5), and we adopted this solution. 

Fig 12 summarizes the results of our de-reddening proce- 
dure, by comparing the derived values of E(B-V) for our com- 
plete sample with the corresponding "optical" values, (B-V)- 
(B - y)o, as a function of Teff (with (B-V) given by the ref- 
erences in Table [Q entry "ref 2", and the intrinsic colors as 
discussed above). 



From this figure, we find no obvious trend of the difference 
in E(B-V) as a function of Teif (the average differences being 
almost exactly zero for supergiants and -0.01 mag for the re- 
maining objects), which is also true if we plot this quantity as a 
function of Mv (not shown). The majority of these differences 
are less than 0.02 mag, which seems to be a reasonable value 
when accounting for the inaccuracy in the observed (B - V) col- 
ors, the uncertainties in the intrinsic ones, the errors resulting 
from our flux calibration and the typical errors on the theoreti- 
cal fluxes (cf. Sect. 13. 2> . 



3. Simulations 

In this section, we will describe our approach to calculating 
the various energy distributions required for our analysis, and 
our approximate treatment of wind clumping, which is based 
upon the assumption of small-scale inhomogeneities. Since this 
treatment consists of a simple manipulation of our homoge- 
neous models, we will start with a description of these. 

Because of the large number of parameters to be varied (M, 
/3, clumping factors), and accounting for the rather large sam- 
ple size, an "exact" treatment by means of NLTE atmospheres 
is (almost) prohibitive. Thus, we follow our previous philos- 
ophy of using approximate methods, whic h are calibrated b y 
means of our available NLTE model grids dPuls et all 120051) . 
to provide reliable results. Note that these grids have been cal- 
culated without the inclusion of X-rays; the influence of X-rays 
on the occupation numbers and IR/radio opacities of hydrogen 
is negligible (e.g.. |Pauldrach et al..200lN) . whilst their effect on 
helium (through their EUV tail) has not been investigated in 
detail. From a comparison of models with and without X-rays 
though, any effect seems to be small. 

We have been able to design interactive procedures (written 
in iDL acting as a wrapper around FORTRAN-programs), which al- 
low for a real-time treatment of the problem, where all required 
fits and manipulation of Hq, spectra and IR/radio fluxes are ob- 
tained in parallel. 

3.1. Ha 

In the present study, synthetic Hq, profiles are calculated as de- 
scribed in Paper 1. This ap proach bases on the approximate 
treatment as introduced by IPuls et al.l (119961) . updated to ac- 
count for line-blanketing effects. Except for the inclusion of 
clumping, no further modifications have been applied; note in 
particular that we have used the same Hq. observations and 
H/He departure coefficients as adopted in Paper 1. 

On the other hand, for most of our sample stars we have 
quoted (and used, within our de-reddening procedure) wind pa- 
rameters from a complete NLTE analysis, which do not rely 
exclusively on Hq, but also on Heii 4686 and other diagnos- 
tics. Furthermore, the observed Hq profiles used here are dif- 
ferent to those in the corresponding sources, because of the 
variability of Hq (cf. Sect. 12. 1> . Thus, we have to check how 
far the values from the complete analysis (denoted by M(in) 
andj6(in)) might deviate from solutions resulting from our sim- 
plified method, used in combination with our different Hq data. 
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Fig. 3. Consistency check for Ho, profiles: results of our approximate H^ line synthesis, for some representative cases from 
Tablels] Dotted: ll„ line profiles with parameters M(in) andy6(in) as derived from a complete NLTE analysis (cf. TableQ; bold: 
corresponding profiles with M - Mi or /3 - P2 (see Table|5}- 



to obtain consistent initial numbers for the following investiga- 
tions and to re-check the reliability of our approach.** To this 
end, we have re-determined mass-loss rates and velocity expo- 
nents, using our observational material, the stellar parameters 
from Tableland the approximate Hq, line synthesis as outlined 
above. Table|5]summarizes the results from this exercise. 

For three objects (CygOB2#7^ HD 14947 and 
HD 210839), no modifications were required at all, whereas 
for the other stars small variations of M were sufficient to re- 
produce our observational data, mostly by keeping the nominal 
velocity exponent. The average ratio between modified and 
input mass-loss rates was 1.07 + 0.22. 

In some cases (particularly for objects with H„ in absorp- 
tion), a second solution is possible, and in all but one case, we 
kept the nominal mass-loss rate constant, while varying /3 (en- 
try ^82 in Table 13 ■ All derived velocity exponents still lie in the 



** Concerning those (four) objects with wind parameters taken from 
Paper I, we have convinced ourselves that the corresponding fits could 
be reproduced. 

' The second solution with /? : 
tion trough. 



0.9 gives a better fit for the absorp- 



expected range. For representative cases, Fig.|5]displays the re- 
sults of our line synthesis, both for models with the nominal 
values, M(in) and j6(in), and for the best-fitting models from 
TablelU with Ml orySa- 

In conclusion, our simplified routine delivers reliable num- 
bers and thus can be used in our further approach to derive 
constraints on the clumping factors. 



3.2. Infrared fluxes 

For the calculation of the infrared fluxes, we closely followe d 
the approximations as outlined by Lamers & Waters (Il984al) . 
with Gaunt factors from .Waters & Lamers. 1.1984.) . The major 
difference concerns the fact that the rad iative transfer is solved 
by means of the "Rybicki algorithm" ( Rvb ickiL 1 197 11) . to ac- 
count for electron scattering in a convenient way. 

A further modification regards the photospheric input 
fluxes which were chosen in such a way as to assure that the 
emergent fluxes, on average, comply with the results from our 
detailed NLTE model grids. 
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Table 5. Consistency check for Hq, mass-loss rates and velocity 
field exponents, for those objects with stellar and wind param- 
eters derived from a complete NLTE analysis (cf. Table 0. All 
mass-loss rates are in units of of lO^^Mg/yr. Mj is the mass- 
loss rate as derived from our approximate method, adopting 
j3\ - y6(in), where possible. In some cases, a second solution 
(Mi,/^!) is possible, mostly for objects with Hq, in absorption 
(see text). 



Star 


M(in) 


yS(in) 


M, 


ySi 


M2 


P2 


CygOB2#7 


10.61 


0.77 


M(in) 


/3(in) 


9.5 


0.90 


HD 15570 
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M(in) 


0.95 
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1.00 
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0.74 


13.00 
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M(in) 


0.95 


HD 210839 


7.95 


1.00 


M(in) 
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6.22 
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5.70 


1.14 


M(in) 


1.05 


HD24912 


2.45 


0.80 


4.00 
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0.80 


1.30 
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M(in) 


0.92 
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1.05 


0.80 


1.30 


/3(m) 


M(in) 


0.90 
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3.07 


1.15 


2.40 
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After some experiments, it turned out that the best choice 
for the various parameters is the foUowing:- 



The velocity law is specified by 

V(r) = V.o(l - b/rf, b^\- (Vmin/Vco)'^^, 



(3) 



where r is calculated in units of /?», and the minimum velocity, 
Vtnin, is set to 10 km s"' . 

Electron temperature. All Gaunt factors are calculated at a 
temperature of 0.9 T^s, and the electron temperature is calcu- 
lated using L ucy's temperature law for spherical atmospheres 
( ILucm.I19711 his Eq. 12, and using grey opacities), with an op- 
tical depth scale accounting for electron scattering only and a 
temperature cut-off at 0.5 Teff . Remember that the radio fluxes 
are almost independent of the temperature, and a number of 
tests have shown that different (reasonable) temperature strat- 
ifications have negligible effects on the derived IR fluxes as 
well. 

Ionization equilibrium. Hydrogen is assumed to be (almost) 
completely ionized, helium as singly ionized outside the re- 
combination radius (see below) and the CNO metals as either 
two or three times ionized. 

Throughout the parameter range considered here, helium is 
singly ionized in the radio emitting region (for A> 2 cm; con- 
cerning mm fluxes see below), as we have convinced ourselves 
by an inspection of our model grid. (Only for 03/4 dwarfs and 
earlier types - which are missing in our sample - does helium 
remain completely ionized throughout the entire wind). 




4.40 4.45 



4.70 



Fig. 4. Location of Hem recombination in velocity space (in 
units of Vco), as a function of Teff. Crosses display this loca- 
tion as derived from our model grid, with 27.5 kK < T^s < 
47.5 kK, diff'erent gravities and wind densities. Curves indi- 
cate the results of our linear regression, Eq.0] for the (limiting) 
cases (log g=3.0,p = -13, solid), (logg=3.0,p = -11, dotted), 
(log §=4.5, p = -13, dashed) and (log §=4.5, p = -ll,dashed- 
dotted). For units of p, see text. The only region which is not 
matched by our regression is the low gravity, low wind density 
region around Teif = 33,000 K, where the regression yields too 
low recombination velocities. 

With respect to the mid- and far-IR emitting region, this 
statement is no longer justified, and one would have to cal- 
culate a consistent ionization structure, which is beyond the 
scope of this paper. In order to obtain an approximate solu- 
tion of this problem, we have parameterized the recombina- 
tion radius, in dependence of Teff, \ogg and mean wind den- 
sity, p = M/{R^Voo), from a linear regression to correspond- 
ing results of our model grid (cf. Fig.|4] crosses). It turned out 
that a best fit could be obtained for the recombination velocity 
(defined as the velocity where the ionization fraction of Hem 
becomes larger than the fraction of Hen when proceeding from 
outside to inside), which can expressed (in units of Voo) as 

V'rec = C -H ai log Teff -H fl2 log g + ai, logp, (4) 

min(vrec) - 0, max(vrec) = 1, 

with 

27,500K<reff <35,000K: 

c = -34.60 ai = 7.79 02 = -0.3325 a^ = -0.0854 

37,500K<reff <47,500K: 

c = -14.90 fli = 3.31 fl2 = -0.0956 ^3 = -0.0798, 

where T^ff is measured in K, M in Mg/yr, R, in Rq and v^o in 
km s '. For 35kK < Teff < 37.5kK, we have applied a linear 
interpolation. 

Concerning the models of our grid, this relation results in a 
mean difference, (yrec(Eq.|3-Vrec(model)) = 0.011 + 0.079 (in 
units of Vco), where the largest discrepancies are found in the 
low gravity, low wind density region around T^ff - 33,000 K 
(cf. Fig. 13 . Note that for high gravity, \ogg = 4.5, and low 
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wind density, logp = -13.0 (dashed line), the complete wind 
contains solely Hen for Teff ^31 ,000 K, whereas for low grav- 
ity, \ogg = 3.0, and high wind density, logp - -11.0 (dot- 
ted line), it remains completely ionized for T^f^ k, 42,000 K. 
For our final model of f Pup (HD66811), our approximation 
yields Viec = 0.87, whic h is in good a greement with the value 
of Viec = 0.83 found bv lHillier et al.l (.1993.) in their paper on 
the X-ray emission of this object. 

Mostly affected by the presence of Hem (compared to the 
assumption that helium is singly ionized throughout the wind) 
is the mid and far IR-band, where the effective photosphere 
might be located below the recombination radius. (In the near- 
IR, the emitted flux is still dominated by the "real" photo- 
sphere.) Except for a few objects, the former wavelength range 
has not been observed so far, so that our predictions remain to 
be verified in the future. Note finally that from the scaling re- 
lations provided by, e.g., Lamers & Waters ( 1984a), the differ- 
ence in the derived mass-loss rates (using Hem instead of Hen 
as the major ion, i.e., no recombination at all) would result in 
a factor of roughly 0.85 for solar helium content. Further com- 
ments on the influence of the helium ionization balance will be 
given in Sect. 14. II 

Photospheric input fluxes were chosen as follows: For A < 
1 fim, we used Kurucz fluxes, whereas for higher wavelengths 
we used Planck functions with riad= 0.87 T^^fi for 1 /im < A < 
2fim, ri-ad= 0.85 Tgff for 2/im < A < 5 yum and ri-ad= 0.9 Teff 
elsewhere. Note that for considerably larger wavelengths, the 
emergent fluxes become independent of the input fluxes, due to 
increasing optical depths. 

We have compared the fluxes resulting from this simplified 
model with those from our NLTE model grid as calculated by 
FASTWiND, for the wavelength bands V to Q. (A comparison be- 
yond 30 fj.m is not possible, since this is the maximum wave- 
length considered in fastwind, which follows from the con- 
straint that, for all IR wavelengths and all wind densities, the 
wind plasma should become optically thick only well inside the 
outermost radius point, /?max = 100 Rf) 

For this comparison, 204 models within the range 30 kK 
< Tgff < 45 kK, with different gravities and wind densities (cor- 
responding to log Q = l og(Q7vi;^)_= -13 .15. . . -12.1, if 



IPuls et alSmOSL 



calculated in km s" , see 
used. As a result, the mean ratio of IR fluxes from our simpli- 
fied model to those from fastwind is of the order of 0.99. . . 1.01 
(different for different wavelengths), and the typical standard 
deviation for each wavelength band is below 5%. 

3.3. Radio fluxes 

Radio fluxes are calculated in analogy to the IR fluxes (with 
identical parameters), but neglecting electron scattering. We 
use a numerical integration, with R^-^x - 10,000 7?*'°, for 
the following reasons: first, the a nalytical express i on by anal- 
ogy to Eq. |2l as provided by IPana^ia & Fellil 1119751) and 



IWright & Barlow! ( 119751) . is valid only under the condition that 
the plasma is already optically thick at v(r) «; Voo, which is not 
the case for objects with thin winds. Secondly, the inclusion 
of depth-dependent clumping factors requires a numerical inte- 
gration anyway. Of course, we have checked that for constant 
clumping factors and large wind densities, the analytical re- 
sults are recovered by our approach. Remember that the emitted 
fluxes are almost independent of the assumed electron tempera- 
ture. From our final results, it turned out that except for the mm 
fluxes of our hottest objects, CygOB2#7 and HD 15570, the 
radio photospheres of the complete sample (even if sometimes 
below v'co) are well above the corresponding recombination ra- 
dius (cf. Table^}. Thus, unless explicitly stated otherwise, he- 
lium is adopted to be singly ionized in our radio simulations. ' ' 
In the following figures, the radio range is indicated to start at 
400 fim - 0.4 mm (end of IR treatment at 200 i-im), but this 
serves only as a guideline, since at these wavelengths helium 
might still not be completely recombined. 

3.4. Inclusion of wind clumping 

To account for the influe nce of wind c lumping, we follow the 
approach as described bv lAbbottet al.l (ll98D. Modified by one 
additional assumption (see below), this approach has been im- 
plemented into NLTE model atmospheres already bv lSchmutJ 
( 1995), and is presently also used by the alternative NLTE code 
CMFGEN. In the following, we will recapitulate the method and 
give some important caveats. 

Regarding the hydrodynamical simulations of radiatively 
dr iven winds, th e term "clumping factor" has been introduced 
bv lOwocki et aLLl988.) . as defined from the temporal averages 
in Eq. ^ To allow for a translation to stationary model atmo- 
spheres, one usually assumes that the wind plasma is made up 
of two components, namely dense clumps and rarefied inter- 
clump material, in analogy to snapshots obtained from the hy- 
drodynamics. The volume filling factor, /, is then defined as 
the fractional volume of the dense gas, and one can define ap- 
propriate spatial av erages for densities and density-squares (cf. 
lAbbottetal.ll98lh . 



<P> 



Sect. 10) have been < P 






f(p^f + (l-f)(p-f]dV, 



(5) 



(6) 



where p^ and p" denote the overdense and rarefied material, 
respectively. Here, and in the following, we have suppressed 
in our notation any spatial dependence, both of these quantities 
and of/. The actual mass-loss rate (still assumed to be spatially 
constant, in analogy to the temporal averaged mass-loss rate 
resulting from hydrodynamics) is then defined from the mean 
density. 



M - A-nr^ 



<p> 



(7) 



'" Within our procedure, we always check that the plasma remains 
optically thin until well inside the outermost grid point. 



and any disturbance of the velocity field (e.g., influencing the 
line-transfer escape probabilities; see .Puis et aL.1993a) is ne- 
glected. 

' ' Concerning the influence of the adopted He ionization on derived 
mass-loss rates, see also;Schmutz & Hamann 1986. 
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The modification introduced by ISchmuta ( 1 1 9951) relates to 
the results from all hydrodynamical simulations collected so 
far, namely that the inter-clump medium becomes almost void 
after the instability is fully grown, i.e, outside a certain radius. 
In this case then, p" — > 0, and we find, assuming sufficiently 
small length scales (see below). 



■■p>-^SWM''-fp^ 



<p'> 



h\WA 



dV = /(p+) 



2 _ <P> 



f 



(8) 



(9) 



Comparing with Eq. [2 and identifying temporal with spatial 
averages, we obtain 



1 

/ci - T 



and 



, < P > 

P = —r- = /ci < P >' 



/ 



(10) 



i.e., the clumping factor describes the overdensity of the 
clumps, if the inter-clump densities are negligible. 

Concerning model atmospheres and (N)LTE treatment, this 
averaging process has the following consequences :- 

- Since, according to our model, matter is present only in- 
side the clumps, the actual (over-)density entering the rate 
equations is p^ - f^x < p > (where the latter quantity is 
defined by Eq. 0. Since both ion and electron densities 
become larger, the recombination rates grow, and the ion- 
ization balance changes. As a simple example, under LTE 
conditions (Saha equation), and for hot stars, we would 
find an increased fraction of neutral hydrogen inside the 
clumps, being larger by a factor of /J compared to an un- 
clumped model of the same mass-loss rate. Further, more 
reahstic, examples for important ions have been given by 
lBouretet"ai](l2005l) . 

- The overall effect of this increase, however, is somewhat 
compensated for by the "holes" in the wind plasma, since 
the radiative transfer (and, consequently, the ionization and 
excitation rates) is affected by the averaging process as 
well, at least for processes which depend non-linearly on 
the density. Note that for processes which are linearly de- 
pendent on the density (e.g., resonance lines of major ions), 
the optical depth is similar in clumped and unclumped 
models, provided that the scales of the clump/inter-clump 
matter are significantly smaller than the domain of integra- 
tion. For p^-dependent processes, on the other hand, the op- 
tical depth is proportional to the integral over < p^ >= /ci < 
p >^^ /ci(p™'^')^, i.e., the optical depths are larger by just 
the clumping factor. Consequently, mass-loss rates derived 
from such diagnostics become lower by the square root of 
this factor, compared to an analysis performed by means of 
unclumped models. 

Before we now comment on the implementation of this pro- 
cess into our models, let us give two important caveats. Implicit 
to the assumption of small length scales, the simple approach 
as described above breaks down (at least to some extent) if 
the clumps become optically thick. In this case, the so-called 
"porosity length" becomes important, and the distribution and 



shape of the clumps has to be specified to allow for more quan- 
titative conclusions. For opacities scaling linearly with density, 
lOwocki et al. (2005) have provided a suitable formalism to de- 
scribe the effects of clumping/porosity in this context, whereas 
for p^-dependent opacities such an analysis is still missing. 

Besides the questions of the length scales involved, related 
optical depth effects and the neglect of velocity disturbances, 
the other important assumption concerns the treatment of the 
inter-clump matter as being void. This approximation is legiti- 
mate as long as clumping is decisive only in those parts of the 
wind which are significantly separated from the base. Under 
this condition, the line-driven instability has already passed its 
linear phase and shocks have developed, compressing the ma- 
terial into clumps and rarefying the medium in between. 

As has been discussed in Sect. ^ recent evidence indicates 
that clumping becomes important from close to the wind base 
on cBouret et al.-.2005.) . In this region, however, the instabil- 
ity is still in its linear phase and resembles more a fluctuation 
(with similar positive and negative density amplitudes) than 
a clumped structure.'^ Consequently, the assumption that the 
inter-clump medium is void becomes questionable. In such a 
case, it might be more appropriate to follow the original ap- 
proach bv Abbott et al. (1981) . namely to account explicitly for 
the "under-dense" medium. 

With respect to our models now, the inclusion of clumping 
effects in the spirit as described above becomes very simple. 
Since all opacities entering the calculations (bound-free, free- 
free and the Hq. line opacity) are dependent on p^, they are mul- 
tiplied with a pre-described clumping factor, whereas the cor- 
responding source functions remain free from such a manipu- 
lation, which is also true for the electron scattering component, 
being proportional to p. Despite our caveats, we assume (/) the 
clumps to be optically thin in H„ and the IR/radio continuum, 
and (//) the inter-clump matter to be void, since, anticipating 
our following results, there is no need to require the inhomo- 
geneity to start already from the wind base on. 

In summary, our procedure is equivalent to other ap- 
proaches used in the literature (e.g., any analysis performed 
with cmfgen), so that the results can be easily compared. 

Since we want to obtain constraints on the radial stratifi- 
cation of the clumping factor, it would be dangerous to use a 
pre-prescribed law, and to adapt only the parameters of such a 
law. An optimum solution would leave the run of the clumping 
factor completely unconstrained, and would derive this quan- 
tity at all depth points from a maximum likelihood method (or 
other optimization algorithms) by fitting to the observed data. 
In view of our interactive procedure, and particularly because 
of our desire to also elaborate on the allowed range of the vari- 
ous possibilities'^, we follow a simplified philosophy, by defin- 
ing five different regions of the stellar wind with corresponding 
average clumping factors, denoted by 



'^ This should be true, even if a different, unknown instability were 
responsible for the development of an inhomogeneous structure. 

'^ Note that, e.g., the velocity-law-index,/?, and the run of the clump- 
ing factor are interrelated, and that for most of our objects observa- 
tional data in the far-IR are missing. 
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The boundaries of these regions and the clumping factors can 
be adapted within our procedure. The first region with fixed 
clumping factor, /d - 1, has been designed to allow for a 
lower, unclumped wind region, in accordance with theoretical 
predictions and our argument from above (namely that any in- 
stability needs some time to grow before significant structure 
is formed). But note also that by choosing r\a = 1 we are alter- 
natively able to simulate a wind where the medium is clumped 
from the wind base on. 

Typical values for rin, r^a, rout and rfar are 1.05, 2, 15 and 
50, respectively. For not too thin winds, this coiTesponds to the 
major formation zones of H„ (region 1 and 2), the mid-/far-lR 
(region 3), the mm range (region 4) and the radio-flux (region 
5). Note that for a number of test cases we have used different 
borders, and sometimes combined region 4 and 5 into one outer 
region. All clumping factors derived in the following are aver- 
age values regarding the different regions, which admittedly are 
rather extended. In almost all cases, however, with such a low 
number of regions consistent fits could be obtained, with rather 
tight constraints on the global behaviour of the clumping factor. 

As a final comment, we like to stress a fact which has been 
mentioned already in Sect. [0 Since (except for electron scat- 
tering) all diagnostics used in this investigation have the same 
dependence on the clumping properties, we are not able to de- 
rive absolute values for the clumping factors, but only relative 
numbers. Note at first that in the case rin = 1 all results de- 
rived for fc\{r) could be multiplied with an arbitrary factor, if 
in parallel the mass-loss rate were reduced by the square root 
of this value, without any loss in fit quality. The only physical 
constraint is the requirement that the minimum value (regard- 
ing all five regions) of the derived clumping factor must not be 
lower than unity. The corresponding mass-loss rate is then the 
largest possible one. 



If, on the other side, r^ + 1, this scaling property is no 
longer exactly preserved, because of the presence of an un- 
clumped region not affected by such a scaling. Since partic- 
ularly the innermost core of H„, but also the optical/near-lR 
fluxes (cf. Sect I2.6l l. are formed in this region, they conse- 
quently deviate from this scaling. As it turned out from the 
analysis performed in the next section, these deviations remain 
fairly small, so that, unfortunately, the derivation of absolute 
values for /d and M will require the use of different diagnos- 
tics. 



4. Constraints on the clumping factor: a combined 
Wa, IR and radio analysis 

4. ^ . T\NO prototypical test cases: ^ Pup and HD 209975 

In this section, we will discuss two prototypical cases in some 
detail before presenting the results for our complete sample. 
We will consider ^ Pup as a representative for a high-density 
wind, with Hq. in emission (this star has the best wavelength 
coverage available within our sample, including fluxes at 25, 
60 yum, 0.85, 1.3 mm and 20 cm), and HD 209975 as a repre- 
sentative for a moderate-density wind (Hq, in absorption). 



f Pup. In the following, we will usually display the results of 
our simulations as done in Fig. |5] namely comparing the ob- 
servations and simulations for H„ in parallel with the IR/radio 
range. Fig. |5] immediately shows the dilemma typical for all 
our objects with H„ in emission: the best fit for H„ requires 
a mass-loss rate typically twice as large as for the radio do- 
main, if homogeneous models are used. The far-IR fluxes are 
also closer to the low-M solution than to the H„-fitting one. 
Let us point out already here that this finding is in agreement 
with a recent comparis on of consistent''* Hp, and radio mass- 
loss rates performed by Fulle rton et alJ (12006 '. who found the 
same factor-of-two discrepancy for a large number of objects. 



i.e., using identical stellar parameters and distances. 
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Table 6. Clumping factors, boundaries of different regions and 
mass-loss rates (in units of lO^^Mo/yr) for three equally well 
fitting models of ^ Pup (strong wind, H„ in emission, (3 - 0.7) 
and for our best fitting models of HD 209975 (moderate wind, 
Hq, in absorption, yS - 0.9). The first solution for ( Pup (which 
optimizes H„) is displayed in the following figures, whereas the 
second one is almost indistinguishable from the first, though 
slightly worse in Hq,, and slightly better in the mid-/far-IR and 
mm range. In the 3rd model it is assumed that helium remains 
doubly ionized everywhere. Note the large difference between 
the clumping properties of the two stars. 
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The derived radio mass-loss rate is considerably larger than 
the corresponding result from Blomme et al. (2003> (using the 
same data set), due to different parameters (larger distance 
and larger helium abundance adopted here). With identical pa- 
rameters, on the other hand, we obtain similar results, M = 
3.7 -lO^^Mo/yr, compared to 3.5 -lO^^Mo/yr. Note also the 
(small) flux excess in the mm-range (with respect to the radio 
fluxes from a smooth model, dotted line), in agreement with the 
findings by Blomme et al. 

Table|6land Fig.|6l(bold), on the other hand, display our best 
solution for a clumped model which consistently reproduces H„ 
and the complete IR/radio band in parallel. In the spirit as out- 
lined above, the mass-loss rate has been chosen from the region 
with lowest clumping, which in this case is the radio domain. 
By setting f^^ to unity then, the adopted mass-loss rate is the 
largest possible one and coiTesponds to the "homogeneous" ra- 
dio mass-loss rate, M= 8.5 ■10"'' Mo/yr, cf. Fig.|5j right panel. 
In this case, the H^j-forming region displays a typical clumping 
factor of 5.5 (from r = 1.12 to 1.5) to 3.1 (from r = 1.5 to 2), 
and yS has been adapted to 0.7 to provide a perfect H^ fit. 

Fig.Qdisplays the change in H„ when a different onset of 
clumping was chosen. If ri,, were 1.3 (dashed profile), the cen- 
tral emission would be missing, whereas for rin = 1 .0 (dotted 
profile, corresponding to a model which is clumped from the 
wind base on), the absorption trough is not perfectly repro- 
duced: the position of maximum depth is located at too high 
velocities, and the trough becomes too broad, resembling our 
best solution for the homogeneous model. 

From our arguments given at the end of Sect. 13.41 it should 
be clear that in particular the latter solution is not unique, since 
an alternative model with all clumping factors multiplied by an 
arbitrary factor /, in parallel with a mass-loss rate reduced by a 



factor of 1/ ^ff, would result in an identical fit. If, on the other 
hand, the perfectly matched absorption trough for our model 
with Tin = 1 . 12 were actually due to a clumping -free lower wind 
base (and not coincidentally matched due to somewhat erro- 
neous departure coefficients and/or the specific observational 
snapshot'^), such a scaling would no longer work (because of 
the presence of an unclumped region), and our solution would 
become "almost" unique, at least regarding the clumping prop- 
erties of the inner wind. 

The "almost" refers to the fact that a different distribution 
of the individual regions, combined with somewhat different 
clumping factors, gives fits of similar quality. The 2nd entry 
of Table|6lis such an example. In this case, we have combined 
the region between r = 1.12 to 2 into one region, whereas we 
have split the outer region, beyond r - 15, into two regions, 
with a border at r = 50. To fit H(j, (with a slightly worse qual- 
ity than displayed in Fig. |6|l, the innermost clumping factors 
had to be reduced (from 5.5 and 3.1 to an average factor of 
5.0), whereas, by adapting the clumping factors in the middle 
and outer part, the fit quality at 60 yum becomes perfect and 
the quality at 0.85/1.3 mm remains preserved Note, however, 
that the overall stratification of the clumping factors is rather 
similar. 

Fig-Eldisplays the advantage of fitting Hq, and the IR/radio 
range in parallel. Although the primary formation region of 
Ha is below 2 R,, it also remains sensitive to variations of 
the clumping factors in the intermediate wind, r $ 15, as can 
be seen from the reaction in the line wings if /d is doubled 
from 2 to 4 (dashed-dotted profile). Of course, a variation of 
the clumping factors in the inner regions (dotted and dashed) 
has even more impact. On the other hand, as displayed in the 
right panel of this figure, the IR/radio band reacts complemen- 
tarily to variations beyond r = 2, although only from the mid-IR 
on (A k 10 /urn). Thus, a combined analysis is able to provide 
tight constraints on the largest possible mass-loss rate and to 
scan the complete stratification of /ci(r) (at least differentially, 
i.e., modulo a constant factor) if the far-IR is well observed. 
Concerning the possible degeneracy of clumping factors and/3, 
we refer the reader to Sect. 14. 31 

FigEl finally, displays the possible error if the helium ion- 
ization were different to that assumed here (cf. Sects. IT2l and 
13.31 ) If helium were singly ionized throughout the complete 
wind (instead of recombining only at v,ec = 0.86), the synthetic 
10 and 20 yum fluxes in particular would become too low; com- 
pensating for this effect by increasing clumping factors is not 
possible, because Hq, would then no longer be fit. If, on the 
other hand, helium were to remain doubly ionized in the out- 
ermost region also, the radio/mm (and the far-IR fluxes) would 
become larger than observed; in this case, a reasonable fit is 
still possible, by lowering the mass-loss rate and increasing the 
inner clumping factors (with a factor roughly corresponding to 
(Moid/Mnesx)^)- The parameters for such a model (which fits 
both Hq, and the entire IR-radio range) is given in Table|6j 3rd 



Co ncerning the temporal variability o f H„ in j ^ Pup, see 



Rfii^^Howarth ( 1996) and references therein. lPuls et alJ il993b) and 
Berghofer et alJ il996'). From these data-sets, a moderate variability of 



the absorption trough is visible indeed. 
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entry. The rather large difference in the resulting (maximum) 
mass-loss rate (factor 0.7) and clumping factors is due to the 
fact that our model of ( Pup has a helium content which is 
twice solar, Ihc = 0.2. For solar helium abundance, as is typi- 
cal for most of the other objects of our sample, the correspond- 
ing factor would be 0.85, as outlined in Sect. 13.31 Note again, 
however, that it is rather improbable that helium is still doubly 
ionized in the radio-forming region. From the consistency of 
the mm and radio fluxes, it is also clear then that the Helium 
ionization must be similar in the mm and the radio forming 
region, in agreement with our predictions for Vrec • 

HD 209975. Table|6landFig.|9ldisplay the results of our com- 
bined fit procedure for this star, which has a moderate wind 



Fig. 8. Clumped models for ( Pup in the IR/radio band: influ- 
ence of helium ionization. Solid: best fitting model, witii Hem 
as the major ion for v < Vrec - 0.86 (5.3 /?»), and Hen as the 
dominant ionization stage outside Vrec; dashed-dotted: Hen as 
the major ion everywhere; dashed: Hem as the major ion in the 
radio emitting domain. 

density and H^ in absorption. Again, we have indicated the re- 
sulting profiles/fluxes when the derived clumping factors are 
varied by a factor of two in specific regions, to check for their 
sensitivity. Most interestingly, this object can befitted with al- 
most constant clumping factors throughout the wind, in stark 
contrast to the above example. Indeed, with slightly different M 
andyS, an almost equally perfect fit is possible with all clumping 
factors being unity. If at all, the (homogenous) radio mass-loss 
rate is somewhat higher than the mass-loss rate derived from 
Ha, SO that in this case /^'" is set to unity. 

Note that a moderate clumping factor of 2 for 1 .5 < r < 2 
is still consistent with the data, and that due to missing far-IR 
information (the indicated data denote upper limits derived by 
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Fig. 9. As Fig.|6l but for HD 209975. The best-fitting model (with all clumping factors at or close to unity) is displayed in bold. 
Other curves show the effects of varying, by a factor of 2, the clumping factors in individual regions alone. Dotted: /ci( 1 .05 . . . 1 .5 
R,) 1^2; dashed: /ci(1.5. . .2 R,) 1^2; dashed-dotted: /ci(2. . .15 R,) 1^2; dashed-dotted-dotted: /ci(> 15 R,) 1.3^2.6. Again, 
Hq, remains sensitive to variations below r = 15 /?* (but see text), whereas the far-IR (not constrained by observations) is mostly 
sensitive to variations in the range 2. . . 15 Rt. Note that the dashed solution is also consistent with the observations. 



IRAS), the clumping in the intermediate wind remains some- 
what unconstrained. After some experimentation, it turned out 
that the data are also consistent with a moderately clumped 
wind (/ci - 10) in the region 10 < r < 50, or a weakly clumped 
wind (/ci - 2) in the region 3 < r < 50 (not quoted in Table|6j. 
Only for the outermost wind (r > 50), do the clumping proper- 
ties have to be similar to the inner wind conditions. 

Since the innermost wind has the lowest clumping, no state- 
ment concerning its onset is possible within our approach. 
Thus, any scaled solution (/d multiplied with /, M reduced by 
1 / ^ff) provides an equally perfect fit and cannot be excluded. 

The second entry for HD 209975 in Table |6l refers to our 
"standard" division of the different regions used for winds with 
Hq. in absorption, namely r^iA = 2 and rout = 10. This scheme 
accounts for the fact that in moderate/low density winds the IR 
and radio emission is formed closer to the star. As can be seen 
from the best fitting clumping factors, the results do (almost) 
not depend on details of the specific borders. 

In summary, the inner and outer wind of this object have 
similar clumping properties, whereas far-IR observations are 
required to constrain the intermediate region. 

4.2. Clumping properties of the complete sample 

Before discussing the results of our analysis for the complete 
sample, let us point out some general findings, and remind the 
reader that the derived clumping factors are independent of any 
uncertainty concerning radius and distance, since all our diag- 
nostics (Hcf/radio/IR) scale in an identical way with respect to 
these quantities, cf. Eq.|2and corresponding discussion. 

The core of Hq, as a tracer of wind clumping below r x 2/?*. 
Our simulations show that the strength of the core of Hq,, 
whether in emission or in absorption, is quite sensitive to the 
value of the clumping factor in the inner part of the wind, and 



thus can be used to determine this parameter out to distances of 
about r w 2. If one relies on the value of y6 as derived by means 
of unclumped models, the corresponding (average) clumping 
factors are very precise, with an accuracy of roughly 10% (but 
see Sect. 14. 3> . Note particularly that clumping factors f^" of or- 
der 2 or larger are still visible for objects with Ha in absorption 
(see Fig.|9j dotted profile). 



Constraints on the clumping factor beyond r «; IR^,. In addi- 
tion to constraining the clumping properties in the lower wind, 
Hq can even serve as an indicator of wind clumping in lay- 
ers beyond r «; 2 (e.g. Fig. |6l left panel). How much be- 
yond? The answer depends, of course, on the specific wind 
density, but some general statements for stars with Hq in emis- 
sion are possible though. Usually, we found that reducing the 
extent of the intermediate clumped region 3 from rout = 15 to 
about rout = 3 has a noticeable effect on the strength of the Hq 
emission wings. The same is true if the boundary of region 2, 
'"mid = 2, is extended to a value of rmid = 3. The effect be- 
comes visible when the outer boundary moves from rout = 15 
to rout = 5 and is insignificant if rout is set to 8 stellar radii 
instead of rout = 15. 

For those objects with Hq in emission and missing far- 
IR/mm information, in Table0 column 14, we have indicated 
the outermost radius, r^^, to which Hq alone can provide in- 
formation on the clumping factor, on the assumption that the 
region, r^^ < r < rf^, is "unclumped" (or, more precisely, 
has the same clumping properties as region 5). In parallel, we 
also quote the corresponding value, Z™"*, which is somewhat 
larger than the original one (for rout = 15), due to the reduced 
width of region 3. Indeed, for almost all objects, r^^ is of the 
order of 5 /?., except for HD 14947 and HD 192639, where Hq 
provides information only out to 3 /?,. Thus, it is safe to con- 
clude that Hq constrains the clumping factor up to distances of 
r = 3. . .5Ri, if in emission. Note, however, that in some cases, 
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significant clumping in region 4 (from rout to rfai) has an effect 
on Ha, which leads to an additional constraint on the clumping 
in this region. 

For objects with H„ in absorption, on the other hand, the 
intermediate region remains much less constrained (Fig.|9] left 
panel), and we will comment below on the corresponding lim- 
its. 

Table summarizes the results of our simultaneous 
Ho-ZIR/radio analysis for the two objects already discussed in 
Sect. 14.11 and for the remaining ones. We have ordered the 
sample according to Ha profile type and spectral type. For al- 
most aU objects, we have used identical boundaries, rin - 1 .05, 
'"mid = 2.0 and rfar = 50, to obtain comparable results. The 
default values for rout coiTespond to 15 (H„ in emission) and 
10 (Hq. in absorption or of intermediate type), but have been 
adapted where necessary. Detailed comments regarding the in- 
dividual objects are given in AppendixlEl where all fits are dis- 
played as well. 

Overall, our simulations show that for stars with Ha in 
emission, a simultaneous fit of the observed radio fluxes and 
the shape and strength of H„, requires clumping factors which 
are always higher in the H„-forming region than in the radio- 
forming one. For stars with H„ in absorption, the situation 
seems to be different: in most cases, the required clumping fac- 
tors are of similar order in the inner and outer regions, as al- 
ready discussed for the case of HD 209975. Note, however, that 
this preliminary impression is dependent on the actual value of 
j3, a problem which will be discussed in our eiTor analysis fur- 
ther below. 

For all objects quoted with a definite mass-loss rate (and 
not only an upper limit), this value represents the largest pos- 
sible value (for given R^), usually derived from adopting an 
outer, unclumped wind with /^^j" - 1 or, for weaker winds, 
^out _ J These mass-loss rates correspond to the "usual" ra- 
dio mass-loss rate. Only for one object, HD 34656, did the 
maximum mass-loss rate have to be derived from Hq,, since 
the radio regime seems to be more strongly clumped than 
the lower wind, at least if the radio emission is purely ther- 
mal. Remember that the radio and Hq. mass-loss rates for 
HD 209975 are consistent to within the error bars. 

Because all our diagnostics depend on p^, different solu- 
tions with lower mass-loss rates and scaled clumping factors 
are consistent with the observational data to a similar accuracy 
as obtained from our fits, except for the innermost cores of Hq 
(particularly if of P Cygni shape), due to our assumption of 
an unclumped innermost region. As already noted, these devi- 
ations remain very small for the derived values of ri,, though. 

For six objects, the maximum mass-loss rate could not be 
uniquely constrained, and the quoted limits correspond to the 
largest value consistent with the data. In five of theses cases 
(denoted by superscripts "a" and "d"), all radio fluxes are up- 
per limits only, and consequently the derived mass-loss rates 
as well. One object (CygOB2#8A) is a confirmed non-thermal 



'* more precisely, where the optical depth t = 1 is reached along the 
radial ray, not to be c onfused with the so-calle d "ef fective radius" lo- 



cated at T '■ 
h984iJ) . 



: 0.24, e.g. lWrJght & Barlow<fl97^ and lLamers & Water j 



emitter ( iBieging et allll989l) . and the adopted maximum mass- 
loss rate relies on the 2 cm which gives the lowest (radio) M 
within the available data set (see Sect. 12. U . 

Of course, all objects with only upper limits for the ra- 
dio flux(es) might be non-thermal emitters, and our interpre- 
tation depends on the assumption that the radio excess is due 
to thermal emission alone. In addition to these objects, three 
more stars (HDi90429A, HD 34656 and HD 37043 (SB2!)) 
have somewhat peculiar radio fluxes, and might also be non- 
thermal emitters. 

Mostly because of these peculiarities (for more details, 
see Appendix |B}, we have given two possible solutions for 
HD 190429 A, HD 37043 and also HD 24912 in Table com- 
prising a minimum and maximum solution with respect to the 
(relative) clumping properties. For HD 37043 and HD 24912, 
the 2nd entries are the more plausible ones (as discussed in the 
appendix), whereas for HD 190429 A both solutions have sim- 
ilar problems (though the difference is not as large as for the 
other two stars). 

Indicated by their name appearing in bold face, the remain- 
ing objects (six with Hq, in emission, one with intermediate type 
and three with Hq in absorption) have well-constrained clump- 
ing properties, i.e., the derived results are robust if/3 is not too 
different from the values derived or adopted here. 

The latter quantity has been specified as follows. For ob- 
jects with Hq in emission and of intermediate type, we have 
used the values from our unclumped analyses (see Tables^and 
13 wherever possible, i.e., if satisfactory fits could be achieved. 
This turned out to be true in almost all cases, with the no- 
table exception of ^ Pup, where our clumped analysis favours 
a much lower value (/6ci - 0.70) than previously found. For 
most objects with Hq in absorption (except Cyg OB2#8C and 
HD 34656, see AppendixiBl. because of missing constraints we 
used the "standard" value (from hydrodynamical models) of yS 
- 0.9, to obtain at least consistent results. Further consequences 
of this uncertainty are discussed in the next section. 

For those stars where Hq is of P-Cygni shape or displays a 
well-refilled absorption trough, conclusive limits could be de- 
rived regarding the maximum value of r^^, i.e., the maximum 
extent of a potentially unclumped region. In all cases, this re- 
gion lies below 1.2 R,. 

In addition to the derived clumping factors which represent 
the best-fitting solution, we also provide maximum values for 
ymid ^jjj yout \yjjj(;jj ^j-g gj-jjj consistent with our data and can 
be restricted further only by additional far-IR and sub-mm ob- 
servations. For all objects with entries "above" Cyg OB2#8C 
in Table0 /™'' could be constrained from the wings of Hq, ei- 
ther for the entirety of region 3 or, if indicated, at least out to 
r^yj. For the other objects, the wind density is too low to induce 
significant reactions in either Hq or the IR when the clump- 
ing properties in region 3 are changed, such that more definite 
statements are not possible. 

4.3. Errors in the derived clumping factors 

In the following, we will concentrate on the errors introduced 
into the derived clumping factors; errors in the mass-loss and 
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Table 7. Clumping properties as derived from our combined H„/IR/radio analysis. Stars are ordered according to Hq, profile type 
("pt") and spectral type. Entries with name in bold are objects with extremely well-constrained clumping parameters. 
Jeff is given in kK, and M^ is the largest possible mass-loss rate, in units of lO^^'Mg/yr. "ratio" gives the ratio of "clumped" 
mass-loss rate to optical results using unclumped models (cf. Table QJ. fi^\ is the velocity field exponent as derived or adopted 
here, v,ec and r^^^ are the velocity (in units of Voo) and radius where He recombines (see Sect. 13. 2L respectively, and r(T2) is the 
radius where the radio continuum becomes optically thick"' at 2 cm (rrec and r(r2) in units of/?,). 

Clumping factors and boundaries are defined as in Sect. 13.41 For all models, region 1 with /d = 1 (not tabulated) extends from 
r = 1 to Tin = 1 .05, except for HD 668 1 1 where rin - 1.12, and rfar (defining the border between region 4 and 5) has been set 
to 50 Rt always. For objects with H,, in emission or of intermediate type, and missing far-IR/mm data, r^^ (with corresponding 
clumping factor) indicates the maximum radius to which H„ alone can provide constraints on the clumping, on the assumption 
that the outer wind is "unclumped" (see text). For objects with Hq, in absorption, /™^ gives the maximum possible clumping 
factor in region 3, which is still consistent with the data. /™'^ is defined similar to /™x, but for region 4. For comments on 
individual objects and corresponding fits, see AppendixlBl 
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"' only upper limits of radio fluxes available; M maximum radio mass-loss rate. 
*' He assumed to be recombined in radio region (see AppendixlBl. 
'■^ upper limit, since non-thermal radio emitter; M from 2 cm flux. 
"" M from H„, since radio fluxes (upper limits only) give larger value. 



modified wind-momentum rates are dominated by errors in the 
angular diameter and radius, but do not affect the major out- 
come of our investigation. 

Let us first mention that during our detailed fits we found 
no systematic problem concerning an underestimation of the 
A^- and Q-band fluxes, so that at least our absolute flux calibra- 
tion seems to be appropriate (see Sect. 12.4.11 . On the contrary, 
for some objects (g-band: HD 15570; A^-band: CygOB2#ll, 
#10 and HD 207198), these fluxes lie above our predictions for 
the best-fitting model. To investigate this point in more detail, 
however, additional fluxes in the mid- and far-IR are required. 

Uncertainties introduced by tine radio continuum. To deter- 
mine the uncertainty in the derived clumping factors due to un- 



certainties in the observations (e.g., intrinsic errors and/or tem- 
poral variability of the observed radio fluxes), we have varied 
/" and f™'^ by identical factors and adapted M accordingly, 
until the observed radio fluxes could no longer be matched. 
From these experiments, it turned out that the clumping factors 
in the regions traced by H(j, (i.e., below r = 3. . .5) are accu- 
rate (on an absolute scale) to within 20 to 50%, whereas the 
ratio of the clumping factors in the various regions remains 
preserved. Remember that the derived clumping factors scale 
inversely with M^ (radio) oc F^-^, i.e., dfd/fd ~ -l-5SFy/Fy. 
Extreme cases regarding this uncertainty in the radio fluxes 
are HD 190429 A, HD 14947 and CygOB2#l 1 (cf. Table|8] 3rd 
column). 
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Table 8. Upper and lower limits for the clumping factors in regions 2 and 3, corresponding to a variation of ySci as indicated ("used" 
refers to the best fitting values tabulated in Table0 For CygOB2#7, HD 15570 and CygOB2#8A we display the solutions for 
the larger values of Z™"^, which fit H^, but somewhat overestimate the 10 fim fluxes, see Appendix 151. No entries are given 
for Cyg OB2#l 1 and HD 34656 due to the very unclear situation encountered for these objects (see Appendix 151. Usually, the 
minimum value of /d refers to the maximum ofySd, and vice versa. For objects with an uncertainty in M being larger than typical, 
column 3 indicates the corresponding range (in units of 10"^Mo/yr). For entries with purely negative AM, the correction refers to 
the maximum value of ySci; in these cases, the outer wind must also be clumped, with values as indicated by /Jj" . For HD 209975, 
the positive correction refers to /Sa =0.7 with ^^|" = 1 , no correction but ^^|" = 1.3 refers to ySd = 0.9, and the negative correction 
and ff' = 3.5 refers to ySd = 1 ■ 1 ■ 
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The degeneracy of /3 and clumping factors in tine inner wind. 
As noted above, the strength of the core of Hq, is highly sensi- 
tive to the value of the clumping factor in the inner part of the 
wind, below r ^ 2/?*. It is also sensitive to the value of the ve- 
locity exponent, /3, and in a similar way: larger values of both 
/3 and clumping factors lead to more emission in the line core, 
giving rise to an unfortunate degeneracy. Note, howeve r, that 
the well-known y6 vs. M degeneracy (e.g.. |Puls et al. 1996) has 
"vanished", since the (maximum) mass-loss rate is determined 
from the radio regime. Except for the weakest winds (which 
cannot be observed in the radio anyway), the radio fluxes re- 
main unaffected by the shape of the velocity field (cf. Table 
column 10). 

This new degeneracy requires an investigation into the 
question of how far any uncertainty in /? will propagate into 
the errors of /d. To this end, we have varied/? and determined 
the appropriate values of /™ and /™'' such that the quality of 
our H„ fit remained preserved. For profiles with H„ in emission 
and of intermediate type, the minimum and maximum values of 
/3 were taken from those solutions which were still compliant 
with the observed profile shape. For objects with H„ in absorp- 
tion, we used reasonable limits, at /? = 0.7 and/? = 1.1, re- 
spectively. Larger values could usually be excluded from the 
profile shape, whereas in certain cases a lower value (though 
being larger than the physical limit, j3 > 0.5) might still be pos- 



sible. This procedure is somewhat similar to our approach to 
resolving the alternative j3 vs. M degeneracy in homogeneous 
winds, when M is derived from H„ alone (cf. Paper I). 

In Table |S] we have summarized the results of our simula- 
tions. As expected, for stars with H„ in emission, the uncer- 
tainty in /3 is not dramatic. This uncertainty leads to an average 
uncertainty in /J" of about + 30%, whereas for objects with Hq, 
in absorption, much larger uncertainties are possible (factors of 
between 2 and 7), if/? were 0.7 instead of 0.9. 

For most of the objects with Hq, in absorption, a larger value 
of /3 (1.1 instead of 0.9) would have some interesting conse- 
quences. Since for these objects the inner clumping factors are 
of the order of unity foi/3 - 0.9, an increase of /3 cannot be 
compensated for by diminished clumping. Consequently, the 
mass-loss rate must be decreased in this case, to reduce the 
wind emission. Table |S| 3rd column, shows that the required 
amount is of the order of 30%. To stiU obtain a consistent fit in 
the radio domain, /^J"" has to become larger than unity, of the 
order of 2. Thus, ;/ low-density winds were to have a veloc- 
ity field exponent larger than the standard one, the differences 
to the objects with emission profiles would become even more 
pronounced: in this case, the outer region would be even more 
clumped than the inner one. Only if y6 were close to its lower 
limit, would the clumping properties of some of the thin winds 
become similar to those of high-density winds. 
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Fig. 10. Clumping factors, /^'" (region 2), for our sample (cf. 
Table|8}, as a function of the distance-invariant quantity, log Q' 
(Q' - M/R];^, with M the largest possible mass-loss rate, in 
units of Mo/yr and R, in units of Rq). Remember that most 
clumping factors refer to outermost clumping factors set to 
unity. Asterisks: objects with H„ in emission; diamonds: ob- 
jects with intermediate Hq, profile type; triangles: objects with 
H„ in absorption. Black colors: objects with definite maximum 
mass-loss rates (corresponding to bold-face entries in Table0. 
Grey colors: objects with upper limits for M and correspond- 
ing lower limits for /^'". Maximum values of /^'" correspond to 
minimum values of jSci, and vice versa for the minimum values. 
The open triangles with solid error bars display the high-M- 
weak-clumping solution for HD 24912 and HD 37043, and the 
open triangles with dashed error bars the alternative low-M- 
strong-clumping solution for these objects. 



Concerning the resulting uncertainties for /™'' (region 3), 
the situation for Hq, emission type objects is similar as for f^". 
The average minima and maxima lie ~ + 20% below and above 
the best-fitting value of /3. For the objects with weaker winds, 
on the other hand, f^"^ still remains unconstrained, and in all 
cases the upper limits as already quoted in Table0remain valid. 

One last comment. Concerning our model(s) for ^ Per 
(HD 24912), we note in Appendix 151 that large values for the 
clumping factors in region 3 (/™'') are required if the small 
emission humps bluewards and redwards of the Hq. absorption 
trough are to be explained by clumping. If, on the other hand, 
/3 were 0.7 for this object, these humps can be created from 
region 2 alone. 



5. Discussion 

5.1. Clumping properties as a function of wind density 

Fig.^Jdisplays the derived clumping factors for region 2 (i.e., 
the first clumped region) as a function of log Q' = log MjR^^, 
i.e., a quantity which is closely related to the mean wind den- 
sity, but is additionally distance invariant. Remember that in the 
present context M is the largest possible mass-loss rate, and that 
most of the derived factors refer to outermost clumping factors 
set to unity. In other words, they have to be regarded as a mea- 



Fig.ll. As Fig. [H but for the ratio /™'V/eP and objects 
with Ha in emission or of intermediate profile type only. The 
star with the lowest ratio (0.3) is ^ Pup. For the three ob- 
jects with a given interval for Z™"* (Cyg OB2#7, HD 15570 and 
Cyg OB2#8 A, see Table 0, we have used the mean value re- 
garding this interval. 



surement of the clumping properties of the inner wind relative 
to the outermost one. Details of the figure are given in the cor- 
responding caption. 

The most important conclusions which can be drawn from 
this figure are the following. For thinner winds with log Q' % 
-7.5 (a regime which is populated by objects with H„ in ab- 
sorption or of intermediate type, but also by the supergiant 
a Cam), the inner wind seems to be clumped by a similar 
degree as the outermost one, at least if we discard the al- 
ternative low-M-strong-clumping solutions for HD24912 and 
HD 37043 (open triangles with dashed error bars). Note that 
if the latter solutions were the actual ones (and we have indi- 
cated that this is rather possible), then both stars are behaving 
completely diff'erent to the other absorption-type stars. 

On the other hand, for stronger winds (almost all stars with 
emission profiles, plus Cyg OB2#8A), the inner wind seems to 
be more strongly clumped than the outermost one, with an av- 
erage ratio of 4.1 + 1.4. Of course, for this class of objects there 
is also the possibility that we encounter moderately (/^'" w 3) 
and stronger {f™ x 5) clumped lower wind regions, or that the 
degree of clumping decreases again towards the largest wind 
densities. However, due to the restricted number of objects, the 
influence of temporal variations (Sect. 12. 1> and the error in- 
troduced by the uncertainty of the continuum flux level, such 
statements cannot be verified at the present time. 

Fig.lllldisplavs the ratio of clumping factors in the interme- 
diate and inner part of the wind, for objects with H„ in emission 
or of intermediate type; for those objects, this ratio could be 
constrained in a rather robust way. In most cases, the clumping 
properties in both regions are either similar, or the (average) 
clumping factors increase moderately from region 2 towards 
region 3, at most by a factor of 2. Let us reiterate, however, 
that region 3 is rather extended (i.e., local values might devi- 
ate from their average ones), and that we cannot derive definite 
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values for radii larger r^^, ^ 5/?*, except for few cases, because 
Hq. becomes insensitive in this region, and strong constraints 
from the IR continuum are missing. Future observations will 
help to clarify this situation. 

For objects with H„ in absorption, at least upper lim- 
its for the clumping factors in region 3, f^^'^, could be de- 
rived (see Tables0and|S). For three well-constrained objects, 
HD 203064, HD 207198 and HD 209975, these upper Umits lie 
between 1.5 and 2, i.e., they might be twice as large as the 
corresponding values for /", but are still rather low. For the 
remaining stars, the maximum values for Z™"* lie in between 4 
and 25, but only for HD 24912 is a large value actually needed, 
if the observed emission humps are to be interpreted in terms 
of clumping and yS were of order 0.9 or larger (see above). 

Concerning the clumping properties in region 4(15 < r < 
50), finally, definite statements are only possible for those 3 
stars observed in the mm region (see below). For the rest, so- 
lutions with /°j"' = 1 are consistent with the observations, but 
larger values (/™x - 2.. . 20, cf. TableQ are possible as well. 
For HD 190429 A, HD 14947, HD 30614 and Cyg OB2#10, H„ 
still reacts to variations of the clumping factor in region 4, and 
/''?"' could be restricted to values from 2 to 4. Since for weaker 

■^ cl 

winds the radio-forming region can extend into region 4, for a 
number of objects with H^ in absorption, /™' is better defined 
than for the rest, particularly for HD 249 1 2 and HD 3686 1 , with 
fT ^ 2. . .3. 

The best-constrained objects within our sample are ^ Pup, 
HD 15570 and HD 210839, due to IRAS (for ( Pup) and mm 
observations. The first of these objects, ( Pup, displays the 
only notable exception concerning the ratio of /""^ and f^", 
namely that region 3 is much less clumped than region 2. In 
other words, maximum clumping must be close to 2 /?», or even 
lower (cf. Fig. Inland Table|6|l. For this star, the derived clump- 
ing factor for region 4 (extending from 15 to 50 /?») is even 
lower than for region 3: at most, /™' $ 1.8. 

For HD 15570, on the other hand, regions 2 and 3 are simi- 
larly clumped, and the derived clumping factor might increase 
even further towards region 4, with Z™* being 5 to 20 times 
larger than the average clumping in the radio-emitting region. 
In the unlikely case that the wind is not recombined at 1.3 mm, 
even /™' = 1 is possible. For this object, the mm measurements 
from SCUBA are extremely valuable, though the rather large error 
bars leave the situation not as clear as desirable. 

For A Cep (HD 210839), finally, the intermediate region is 
more heavily clumped than the inner one, whereas region 4 
could be constrained (again via scuba observations) to display 
clumping factors between 1 and 8. It remains to be clarified 
whether the two different observed flux levels (Tab.|3 Fig lB.2l 
2nd row) are a sign of significant temporal variability of the 
outer wind (indicating a temporal variation of clumping or a 
non-negligible effect of macro-structures) or the "truth" lies in 
between both measurements, which are still consistent within 
the claimed error bars. 

In summary, at least one of these three objects is rather 
weakly clumped in region 4. Although the same might be 
true for the other two stars (accounting for the lowest possible 
fluxes), a significantly clumped outer region is more probable. 



5.2. Wind-momentum-luminosity relation 

Before discussing some further implications of our findings, let 
us consider the wind-momentum-luminosity relation for our 
stellar sample, accounting for the results derived in the present 
paper. Fig.[21displays two such relations, in comparison with 
the theoretical predictions by Vinket al. (2000).'^ In the left 
panel, we show the results using H„ mass-loss rates derived 
by unclumped models, updated for a re-determined stellar ra- 
dius (Sect.|2j- As already noted in Paper I and outlined in the 
introduction, objects with Hq, in absorption and of intermedi- 
ate type are perfectly consistent with the predictions (except 
for a few objects at log L/Lq < 5.35; see below), whereas ob- 
jects with H„ in emission populate a strip parallel to, but above, 
the predictions. Only the large-distance solution for ( Pup lies 
on the relation, whereas the low-distance solution displays the 
same discrepancy as the other stars (both solutions indicated by 
"ZP"). 

In the right panel, we display our new results, with mass- 
loss rates from Table These mass-loss rates are the largest 
possible ones, and are essentially the radio mass-loss rates if 
the winds were unclumped in the radio-forming regime. Except 
for this assumption, the largest eiTors present in this figure are 
due to errors in the distance estimate. We have deferred from a 
rigorous error analysis concerning this problem, as this is be- 
yond the scope of the present investigation. 

What is obvious from this plot, however, is that the agree- 
ment between observations and theoretical predictions has sig- 
nificantly improved. Almost all objects now lie very close to 
the theoretical relation, independent of profile type. 

The reason, of course, is that the newly derived (radio) 
mass-loss rates for emission-profile objects are smaller than the 
H„ mass-loss rates (see Table0 column "ratio"), by an average 
factor of 0.49 + 0.10. Most interestingly, this is almost exactly 
the same factor which has been claimed in Paper I (0.48, drawn 
from a much larger sample), and which has been used a pri- 
ori in our de-reddening procedure (see Sect. 12. 6>. A factor of 
the same order (0.42) has also been found by iFullerton et alJ 
(2006), for a sample comprising objects similar to those con- 
sidered here. For objects with Hq, in absorption and of interme- 
diate type, Hq, and radio mass-loss rates agree well, and they 
remain at their "old" position. Note that for the only absorption- 
type object in the sample of Fullerton et al. with Hq and radio 
data available in parallel (HD 149757), a comparable agree- 
ment was found, supporting our results. 

Due to the shift in wind-momentum rate, the new position 
of ^ Pup (larger distance) is completely inconsistent with the 
rest, whereas the conventional, lower-distance solution matches 
the predictions perfectly. The same problem was found in 
Paper I (after applying an average down-scaling of wind mo- 
menta, in anticipation of clumping effects), and our present re- 
sult (which confirms this expectation) seems to favour a lower 
radius. 

In accordance with our reasoning in Paper I, however, 
we like to point out that ( Pup is a "bona fide" runaway 



" which are consistent with independent investigations by our 
group, see iPuls et al. (2003) . and also Kudritzki (2002) . 
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Fig. 12. Wind-momentum-luminosity relation for our sample. Modified wind-momentum rate, D^ 
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units. Left panel: mass-loss rates derived from Hq,, using homogeneous models, cf. Table[2 Right panel: largest possible mass- 
loss rates, from this investigation. Upper limits indicate those cases where radio fluxes are upper limits and/or non-thermal 
emission cannot be excluded. Asterisks: objects with Hq. in emission; diamonds : objects wi t h inter mediate profile type; triangles: 
objects with Wa in absorption. Dashed line indicates theoretical prediction bv Vink et alJ (l200r)l) . "ZP" indicates the large and 
low distance solution for [, Pup (see text). For the three objects at log LJLq = 5.3 (HD 203064, HD 207 198 and HD 209975), the 
lowermost solution indicates their position if the velocity exponent was larger than expected (/? = 1.1 instead of /? = 0.9). In 
this case, the (unclumped) Hq. mass-loss rate would be lower than the radio mass-loss rate, and the wind would have to be more 
strongly clumped in the radio regime than in the innermost region. For HD 37043, at log L/L© — 5.45, the lowermost solution 
corresponds to the 2nd entry in Table 



star, (i.e., its parent association. Vela R2, has been iden- 
tified by Sahu & Blaa uw 1993). Based on Hipparcos data, 
IVanbeveren et alJ ( ll99S) argued that f Pup could have become 
a runaway as a result of a supernova explosion in a massive 
close binary, which might explain its peculiar characteristics, 
such as enhanced He and N abundances at the stellar surface, 
high peculiar and rotational velocities, and its overluminosity. 
The reason why the wind-momentum rate should be lower than 
for other objects remains to be clarified though. 

Whereas the "new" WLR agrees extremely well with the 
theoretical predictions for objects with log L/L© > 5.35, the 
three best-defined absorption-type stars at the lower luminosity 
end of our sample (HD 203064, HD 207198 and HD 209975) 
lie too high, by a factor of ^ 2.5. A similar efifect was found 
in Paper I, though at the time it was not clear whether or not 
their winds are clumped in the H„-forming region to a simi- 
lar degree as emission type objects. In addition to wind mo- 
menta based on nominal radio mass-loss rates, we have also 
indicated (by the lower end of the displayed bars) their posi- 
tion if the velocity exponent were to be larger than expected 
(y6 = 1.1 instead ofyS = 0.9). In this case, the (unclumped) H„ 
mass-loss rate would be lower than the radio mass-loss rate, 
and the wind would have to be more strongly clumped in the 
radio regime than in the innermost region (cf. Table |S}. Even 
in this case, a discrepancy of a factor of w 1.7 would still be 
present. To unify these objects with the others by clumping ar- 
guments alone would require that they have to be much more 
clumped in the radio regime (on an absolute level). 

Of course, one might argue that this problem is not related 
to (unknown) physics but to wrong distances and radii. Though 



this might be possible accounting for the mean errors in modi- 
fied wind-momentum rate (0.13 dex) and luminosity (0.19 dex) 
derived for Galactic objects in Paper I, it is more plausible to 
invoke physical reasons, since we have to explain an identical 
problem for three different stars (with different Teff ) at identical 
positions in the diagram. 

Again, we stress that all displayed positions rely on the 
derived, largest possible mass-loss rates. If the radio regime 
were clumped, downward corrections become necessary. In 
this case, however, the displayed agreement would be pure co- 
incidence. 



5.3. Implications and conclusions 

The results from the previous section have confirmed our ear- 
lier hypothesis that the "old" Hq, mass-loss rates for objects 
with Hq in emission are contaminated by clumping, and that, 
compared to theoretical predictions, these mass-loss rates are 
overestimated by a factor of at least 2. . . 3. Regarding the dif- 
ferent behaviour of objects with emission and those with ab- 
sorption profiles, however, we seem to have invoked a some- 
what erroneous explanation (see our arguments recapitulated 
in Sect. H). 

Indeed, if objects with Hq in absorption were clumped in 
the lower-wind region, in a similar way to emission-type ob- 
jects, we would have seen this: note that the presence of clump- 
ing with factors as low as f^" = 2 is clearly visible (cf. Fig.|9}. 
It must be stressed, however, that our present sample consists 
of supergiants and giants only, and that dwarfs (with a very 
low wind emission inside the core of Hq) are missing. At least 
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for the latter luminosity class, our old arguments might still be 
valid. For example, if the wind base was actually unclumped 
(a s allowed for by ou r analysis, but in contrast to the findings 
bv lBouret et al.l2005l) . and H^, predominantly forms in this re- 
gion, we would not see the clumping effects, though they would 
be present in, e.g., the mid-/far-IR. 

Assuming for the moment that, on an absolute scale, the 
outer clumping properties are independent of wind density, our 
results imply that the different degree of consistency between 
the theoretical and observed WLR is likely related to a physical 
effect: inside the H^-forming region, denser winds are more 
strongly clumped than (most of) the weaker winds, at least if/? 
is of the order of 0.9 or larger in the weaker winds (see Fig.llO>. 

What might be the origin of such a difference? Objects with 
H„ in emission have a large wind density and are usually super- 
giants with low gravity and a considerable Eddington-F. Thus, 
it is rather possible that they are subject to photospheric insta- 
bilities and/or pulsations, triggering a somewhat larger struc- 
ture formation in the lower wind, compared to lower density 
winds from higher gravity objects. Indications of such a depen- 
dence are consistent with investigations regarding photospheric 
line-profile variability (increasing with stellar radius and lumi- 
nosity), as outlined by Fullerton et al. ( 1996). 

On the other hand, our findings are in some contrast to hy- 
drodynamical simulations, at least regarding self-excited struc- 
ture formation. If there was any dependence on wind density 
predicted at all, thin winds should be more strongly clumped 
than thick winds, because of the missing stabilization due to 
the continuum (Owocki & Puis, 1999), which induces a more 
heavily structured wind in the lower part. Note also, that in thin 
winds, (transonic) velocity curvature terms become important, 
leading to gradient terms in the source functio ns and modifi- 
cation o f the line acceleration ( Puis et allll99 8). Simulations 
bV iOwocki & Pulsl ([1999i) accounting for this effect resulted in 
a highly structured wind, with a moderately reduced mass-loss 
rate and a rather steep velocity law in the lower region. Thus, 
even the possibility that fi is low (which would increase the 
derived, lower clumping factors, cf. Table |8j cannot be com- 
pletely excluded, although in this case, /™ S 2 (for the three 
well-defined objects) is still rather low. 

Let us now compar e our results with the predicted ra dial 
stratification of /d itself dRunacres & OwockiLl2002ll2005]) . As 
is true for our results, theory also depends on a number of as- 
sumptions. Most important in this context are:- 

- the dimensionality of the hydrodynamical treatment, which 
is mostly 1-D (because of the complexity of calculating 
the radiative line force). First results from a 2-D approach 
(iDessart & Owocki, 2003) might indicate somewhat lower 
(factor of 2) clumping factors than those resulting from a 
1-D treatment. 

- the excitation of the line-driven instability. Almost all mod- 
els investigated with respect to the clumping factor re- 
fer to self-excited perturbations. Unfortunately, externally 
triggered perturbations, suc h as sound waves and photo- 
spheric turbulence (see Fel dmeier et alJll997h . and photo- 
spheric pulsations, have not been examined with regard to 
this quantity. 



Table 9. Clumping factors as predicted by hydrodynamical 
simulations from Runacres & Owocki (2002), for the differ- 
ent regions as used in this investigation. The first average is 
a straight one, the 2nd is weighted with p^ (see text). Note that 
these numbers are only approximate ones, since they have been 
derived from figures and not from tables. 
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- the so-called line-strength cut-off, A-max. In order to keep 
the problem nur nericallv treatable (i.e ., to avoid too fine a 
grid resolutiont. lOwocki et alJ (Il988l) introduced an opac- 
ity cut-off regarding the driving lines, which is typically 
thr ee dex below the actual val ue. Experiments performed 
by iRunacres & Owockil ll2002h showed that the clumping 
factor in the outer wind (around 50 Rt), in particular, can 
increase if more realistic values are used. The inner and 
outermost part seem to remain rather insensitive, at least if 
very low values for /Cmax are avoided. 

Thus, the numbers which will be quoted in the following 
might be considered in a qualitative sense, especially since, 
for our comparison, we have to estimate appropriate spatial 
averages over the different regions. In our approach, we have 
used clumping factors assumed to be spatially constant within 
certain regions, whereas Runacres & Owocki (2002) display 
the clumping factor as a function of r. The most decisive 
quantity regarding radiative transfer is the optical depth, be- 
ing proportional to the spatial integral over fci(r)p(r)^ (assum- 
ing the source function to be unaffected by clumping), so that 
a meaningful comparison requires the predicted clumping fac- 
tors, /ci(r), to be averaged over p^ inside the regions consid- 
ered.''^ To this end, we have used the results displayed in the 
various figures provided by Runacres & Owocki. 

Table 121 summarizes the predictions. Region 1 (the inner, 
unclumped region) typically extends to 1.3 /?, (for thin winds, 
it might be narrower; see above), which is fairly consistent with 
the derived maximum extent of such a potentially unclumped 
domain (rin ^ l.l.. .1.2 Rt). 

Regarding the other regions, we have to discriminate be- 
tween absolute numbers and numbers referring to the average 
clumping factor in region 5, which is of the order of 4. . . 5 or 
even larger, if /Cmax is increased beyond its "standard" value. 
Such large averages depend on results indicating that the outer- 
most wind (beyond 1000 R ,) is also considerably clumped; cf. 
iRunacres & Owockil 112005 '). Only for rather low values of K^nax 
is a smooth radio regime predicted. 



'* By adopting this approach, we discard certain details, such as the 
fact that H„ reacts to averages over constant velocity surfaces (and not 
along the radial direction), as well as optical depth effects. 
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For region 2, we find average values f^" ^ 2. . .3 (lower 
than in region 5!), for region 3 values around 4. . .5, and for 
region 4 values around 11, which again might be even larger 
for large ATmax- Note that for different wind densities and wave- 
lengths, the calculated averages for regions 3 and 4 might be 
higher and lower, respectively, than the indicated ones, depend- 
ing on the radial position at which r = 1 is reached. Finally, the 
predicted maximum is located at the border between regions 3 
and 4 (around 15 /?,), but might be shifted towards larger radii 
for larger /Cmax- 

Compared to our results, these predictions are significantly 
different, at least ;/ the average clumping factor in the radio 
domain is of the order of 4 or larger. In this case, all Hq, mass- 
loss rates should be lower than the radio mass-loss rates, which 
is definitely not true. Thus, either the clumping factors in region 
2 are predicted as too low, or those in region 5 as too large! 

Disregarding this problem, the average clumping factor 
should increase monotonically from region 2 to 4 according 
to theory, and at least some of our emission type objects (e.g., 
HD 15570) are compatible with this result (though for others 
/™' is of the same order or even lower than Z™''). Only con- 
cerning the differential behaviour of region 2 to region 3, do 
most objects behave as predicted (Fig. II U . As outlined akeady 
above, the notable exception to this rule is ( Pup, where the 
complete run of fdir) and the position of its maximum defi- 
nitely deviate from the predictions (and from the ot her objects 
inves tigated). Such a deviation was already found bv lPuls et alJ 
lll99 3b). who tried to simulate the observed H,, profile and IR 
continuum for f Pup, based on hydrodynamical models from 
S. Owocki. Though they were quite successful in fitting H„ 
with a mass-loss rate just a factor of 2 lower than when using 
homogeneous models (and consistent with present estimates), 
the IR continuum was too strong at this M, indicating lower 
clumping factors than predicted in region 3. 

The real question, of course, concerns the absolute value 
of the clumping factors, and their dependence on stellar pa- 
rameters. What has been derived in this investigation is the be- 
haviour of the inner clumping properties relative to the outer- 
most ones. To reiterate, if the outer clumping properties were 
independent of wind density and/or stellar parameters, thinner 
winds would be less clumped in the inner region than stronger 
winds, and we have indicated above a possible reason for this. 
If, on the other hand, the (absolute) clumping factors in the in- 
ner part were to be equal or even larger in thinner winds than 
in denser ones, we would meet a number of other problems 
requiring explanation. In this case:- 

- the outer region of thinner winds has to be more clumped 
than in thicker winds. 

- the consistency with the theoretical WLR would com- 
pletely vanish. 

- the WLR would again show a strong dependence on lu- 
minosity class and/or H^ profile type (even if the theoreti- 
cally predicted off-set was wrong). Such a dependence is 
presently not understandable, since the major prediction 
of radiation-driven wind theory is that the modified wind- 
momentum rate should be dependent on luminosity alone 
(at least if the slope of the corresponding line-strength dis- 



tribution function is not too different from its presently de- 
rived value). 

The only way to clarify this situation is the inclusion of pro- 
cesses which do not depend on p^. One such dia gnostic is Pv 
('Massa et al.', '2003'; 'Full erton et all l2004 12006') which under 
favourable circumstances scales oc p alone. The major problem 
here arises from the uncertainties regarding the ionization frac- 
tion of this ion, which might be additionally contaminated by 
the UV-tail of the X-ray emi ssion. Assuming that P v is a ma- 
jor ion between 04 and 07, iFullerton et alJ 1120061) derived a 
median reduction in M (compared to homogenous Hq, and ra- 
dio diagnostics) by a factor of 20, where thin winds seemed to 
be more affected than thicker ones. Note that this would imply 
clumping factors of the order of 100 in the radio regime! 

Detailed NLTE investigations accounting for clumping, on 
the other hand, are only in their infancy, and again, the inclu- 
sion of X-ray effects is a difficult task. The only object within 
our sample which can be comp ared with suc h an investiga- 
tion is HD 190429A, analyzed bv lBouretet ai] (Eo05). In their 
conclusions, they quote a reduction of a factor of three in M, 
compared to a homogeneous mass-loss rate of 6 -lO^^Mo/yr 
derived from the far-UV, exploiting p- and p^-dependent pro- 
cesses in parallel, and accounting for a consistent ionization 
equilibrium. 

The derived homogeneous UV mass-loss rate is much 
lower than our homogeneous H„ value (radius and distance 
are comparable), an d they specula t e on strong variations in 
Ha, referring also to IScuderi et al.l (119981) . who report an in- 
crease of the H„ equivalent width between 1988 and 1991, by 
a factor of 2 (but see also Markovaet al. 2005. who found no 
indications of such large changes in Hq., at least over an in- 
terval of one year between 1997 and 1998). Though the im- 
plied clumping factor (from a comparison of homogeneous and 
clumped UV mass-loss rates) would be not too different from 
"our" value, on an absolute scale there are much larger differ- 
ences. Comparing their final mass-loss rate (1.8 -lO^^Mo/yr, 
with Rt - 19.5 Rq and Vo, = 2300km s"') with our radio mass- 
loss rate (7.5... 9.5 ■lO"''M0/yr, with R, = 22.7 Rq and v^c 
- 2400km s"'), this would suggest a strongly clumped radio 
regime, with /^^j"" «j 10. . .16, at least if there have been no ma- 
jor changes in the average wind properties between their UV 
and our radio observations. Additionall v. Bouret et al. (2C 
point to the fact that the predictions by iLenorzer et alJ «2C. 
concerning BrQ, indicate that the outer winds "would be less 
affected by clumping", compared to the regions they could ac- 
cess. Thus far, the situation remains unclear. 



Notably, the other object investigated by Bouret et alJ 
(J 2005 ) is an object with Hq, in absorption, and for this object 
they find a reduction in M by a factor of 7 (again with respect 
to UV observations alone). This result would agree with our 
statement from above that thin winds are expected to be more 
strongly structured than thick winds, at least if the latter are not 
externally triggered by photospheric disturbances. 

Accounting for th ese findings and other investigations with 
similar results (e.g., iHillier et all l2003t iBouret et all l2003l) . 
there seems to be increasing evidence that the agreement be- 
tween the theoretical and observed WLR (which, if real, would 
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imply a smooth wind in the radio regime) is indeed just coin- 
cidence, and that the radio regime must be strongly clumped, 
maybe even more strongly than presently described by hydro- 
dynamics. 

Aside from the major implications such a reduction of 
mass-loss rates would have, e.g., regarding stellar evolution in 
the upper HRD and feedback from massive stars, such a result 
would also lead to the following problem: since the present the- 
oretical WLR originates from consistent calculations of the ra- 
diative line force, lower wind momenta would imply that too 
much radiative pressure is available. A reduction of this quan- 
tity, however, is rather difficult (but see below). 

Finally, let us note that a significant down-scaling of 
mass-loss rates would unfortunately also aff^ect stellar pa- 
rameters (again!). For the p^-dependent results derived here, 
such scaling is easily possible, without modifying any result. 
Photospheric lines, on the other hand, might be differently af- 
fected by a strongly clumped, but weaker wind, since they do 
not always scale with Q, but depend on other combinations of 
M, Rt and Vco as well. 

6. Summary and future work 

In this investigation, we have performed a simultaneous anal- 
ysis of H„, IR, mm (if present) and radio data to constrain 
the radial stratification of the clumping factor in a sample of 
19 O-type supergiants/giants, with dense and moderate winds 
(H„ in emission and absorption). All analysis tools used in- 
volve certain approximations, but we have ensured that the 
derived results comply with state-of-the art NLTE model at- 
mospheres, by comparing and calibrating to a large grid of 
such models. Clumping has been included in the conventional 
approach, by manipulating all p^-dependent opacities and as- 
suming the inter-clump matter to be void. Caveats have been 
given to this assumption and other problems inherent to this ap- 
proach, namely the neglect of disturbances of the velocity field 
due to the clumps, and the assumption of small length scales, 
related to the problem of porosity. 

Instead of adapting the dumping-factor at each radial grid 
point (which is possible only if using optimization methods, 
requiring a well-sampled observed wavelength grid), we have 
introduced 5 different regions, with constant clumping factors 
inside each region. Because all our diagnostics depends on p- 
(except for the small contribution by electron scattering), the 
most severe restriction within our approach is given by the fact 
that we cannot derive absolute clumping factors, but only fac- 
tors normalized to a certain minimum. Since in all but one case 
(HD 34656) this minimum was found to be located in region 5 
(or, in other words, since in all those cases the radio mass-loss 
rate is the lowest), our normaUzation refers to the radio regime, 
and the corresponding (radio) mass-loss rate as derived here 
is the largest possible one. Other solutions are possible as well, 
with all clumping factors multiplied by a constant factor, /, and 
a mass-loss rate reduced by ^ff. 

Our analysis is based on H„ line profiles, near-/mid-/far- 
IR fluxes taken from our own observations and the literature 
(de-reddened as detailed in Sect. 12. 6> . mm fluxes observed by 
scuba/sest (own and literature data), and radio data taken from 



our own VLA observations and the literature. We have dis- 
cussed the issue of non-simultaneous observations: based on 
present-day observational facts, the Hq,, IR and radio variabil- 
ity of thermal emitters is low enough so as not to pose any 
problems for our study, at least if the derived results are consid- 
ered in a statistical sense. Within our sample, there is only one 
confirmed non-thermal emitter (Cyg OB2#8 A), and three more 
objects display somewhat peculiar radio fluxes (HD 190429 A, 
HD 34656, see above, and HD 37043). These objects might be 
non-thermal emitters as well, but this has to be confirmed by 
future observations. In any case, the derived mass-loss rates 
(from the minimum radio flux) can be considered as an upper 
limit. 

As it turns out, the core of Hq. provides very useful diagnos- 
tics for the clumping properties in the inner wind (r S 2Ri,), 
and, if in emission, the wings can be used to constrain the 
clumping inside the first five stellar radii, with an additional 
check provided by IR data. If mm fluxes were available, the 
outer wind (15/?* ^ r ^ 507?*) could be constrained as well. 
Only the region between 57?* ^ r ^ 15/?* remains "terra incog- 
nita" in most cases, due to missing far-IR fluxes. 

For ten stars in our sample (six with H„ in emission, one of 
intermediate type and three with Hq, in absorption), the derived 
clumping factors are robust and lie within well-constrained er- 
ror bars. For six stars (including CygOB2#8A), only upper 
limits for the radio mass-loss rate are available, and the de- 
rived clumping factors have to be considered as lower limits. 
Obvious differences to the best-constrained objects were not 
found though, except for HD 24912, which behaves atypically. 
The three remaining objects constitute HD 34656, which is the 
only object in our sample with an H„ mass-loss rate lower than 
the radio mass-loss rate (and as such has been discarded from 
our further analysis), HD 37043, which exhibits similar prob- 
lems to HD24912 (but has a better-constrained radio mass- 
loss rate), and HD 190429 A, which displays a certain degree 
of radio-variability. Taking the various results together, we can 
summarize our findings as follows:- 

- for almost all objects (except for 3 stars with H,, in absorp- 
tion and logL < 5.35 Lq), the derived (radio) mass-loss 
rates are in very good agreemen t with the predict ed wind- 
momentum-luminosity relation JVinket al.. 2000), in con- 
trast to previous results relying on unclumped Hq, data 
alone. If ( Pup is located at the "close" distance, then it be- 
haves as the rest. If, on the other hand, it is located further 
away, its (radio) wind-momentum rate would lie consider- 
ably below the predictions. 

- the mean ratio of radio mass-loss rates to unclumped He 
mass-loss rates for stars with H(j, in emission is 0.49 + 0. 10. 
This is almost exactly the same factor as found in Paper I, 
by shifting the observed WLR (using unclumped models) 
for these objects onto the predicted one. It also agrees well 
with recent findings from lFullerton et alJ OOOfl) . 

- the average, normalized clumping factor in the innermost 
region (r ^ 2/?*) of stars with H„ in emission is ~ 4.1 + 1.4. 

- thinner winds with H^, in absorption have lower normalized 
clumping factors in this region. For all three stars with ro- 
bust constraints, these factors are similar to those in the ra- 
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dio region, at least if the velocity exponent is not too differ- 
ent from the hydrodynamical prediction, /3 x 0.9. Factors 
of the order of P" ^ 2 can be excluded, due to the sensitive 
reaction of H„. 

- for all objects where H„ is of P Cygni shape, or displays 
a well-refilled absorption trough, the maximum extent of a 
potentially unclumped region can be limited to lie inside 

- in most cases, the clumping factors in the inner and adjacent 
region (2/?* $ r ^ 5. . .15/?*) are comparable or increase 
moderately from inside to outside. Only for ^ Pup, does 
our analysis restrict the maximum clumping at r ^ 2/?*. 

- the presence of clumping introduces a new degeneracy 
in the results, namely between the velocity field expo- 
nent, /3, and the clumping factors. If /3 is lower than as- 
sumed or derived from the fits, the clumping factors are 
larger, and vice versa. Extreme deviations of /3 from val- 
ues obtained from an unclumped analysis can be excluded 
though. Interestingly, a perfect fit for ^ Pup requires /3 = 
0.7, contrasted with j3 = 0.9 from unclumped diagnostics 
fRepolustetal.,2004). 

- two of the three stars with mm-observations (HD 15570 and 
HD 210839) indicate a certain probability that the outer re- 
gion 4 (15/?* S r 5 50/?*) is considerably more clumped 
than the radio domain (but remember the rather large eiTor 
bars on the mm data), whereas the third star, ( Pup (with 
negligible observational errors), displays similar clumping 
properties in both regions. 

- Our results differ from hydrodynamical predictions (incor- 
porating the in trinsic, self-excited line-driven instability, 
iRunacres & OwockL200Z.2005.) at least in one respect: the 
latter imply a larger radio than H„ mass-loss rate (or, alter- 
natively, lower clumping in the inner than the outer wind), 
which is definitely not true for our sample. 

In addition to the conclusion that one of the best-observed mas- 
sive stars, ^ Pup, might be a rather atypical representative of 
its kind (maybe due to its possible expulsion from a close bi- 
nary system), the major implications of these findings can be 
stated within three different assumptions concerning the clump- 
ing properties of the outermost regions :- 

assump. (a): The radio region is not, or only weakly, clumped. 
In this case, our "old" hypothesis (concerning a shift of 
mass-loss rates for objects with Hq. in emission, due to 
clumping) would be confirmed, but there would be a phys- 
ical difference between denser and thinner winds, in the 
sense that thinner winds would be less clumped than thicker 
winds in the inner region. This difference might then be 
related to different excitation mechanisms of structure for- 
mation. If assumption (a) were true, the theoretical WLR 
would be perfectly matched. On the other hand, the ab- 
solute numbers for clumping factors and mass-loss rates 
would be in severe contrast to results from other investi- 
gations that have used alternative diagnostics, not directly 
affected by clumping (e.g., the Pv resonance lines). 

assump. (b): The radio region is strongly clumped, but the out- 
ermost clumping factors are independent of wind density. 
In this case, a unification with results from other diagnos- 



tics is possible, and the present mass-loss rates would have 
to be significantly revised, with serious implications for 
the evolution of, and feedback from, massive stars. Again, 
weaker winds would be less clumped in the inner region, 
and the theoretical WLR would no longer be matched. 
One of the most robust predictions from radiation-driven 
wind theory, namely that the modified wind-momentum 
rate should depend almost exclusively on luminosity (and 
not on mass or gravity), would still be consistent with our 
data, even if there were an offset between the theoretical 
and observed WLR. 
assump. (c): The radio region is strongly clumped, but the de- 
gree of clumping is different for different wind densities. 
This case is also consistent with present data, but would 
again imply, in addition to different offsets between the the- 
oretical and observed WLR, that the observed WLR is de- 
pendent on a second parameter. 

Obviously, the implications of all three assumptions pose their 
individual problems, and would have different consequences 
regarding the urgent question about the "true" mass-loss rates 
of massive stars. Since there is no direct way to measure the 
clumping in the radio regime, for further progress we suggest 
the following steps. 

On the observational side, we have to: (/) re-observe some 
problematic objects at radio frequencies, to check their vari- 
ability and to obtain further clues as to whether their emission 
is of thermal or non-thermal origin; and (//) most importantly, 
accumulate far-IR and mm observations, to constrain the (nor- 
malized) clumping factor in the intermediate wind. 

Once a reliable, normalized stratification has been ob- 
tained, it can be used as an input into state-of-the-art model 
atmosphere codes allowing for the inclusion of clumping and 
X-ray emission, with the mass-loss rate/velocity field adapted 
until all diagnostics (including the FUV/UV) are reproduced. 
This would also clarify the question concerning the ionization 
fraction of Pv. After having analyzed a significant number of 
objects, covering a large parameter space, we should be able to 
determine the importance of clumping, how it varies with spec- 
tral type and wind density, and what the actual mass-loss rates 
are. 

Additionally, the derived clumping factors have to be in- 
corporated into stationary wind-dynamics models. Using such 
models, we can investigate how far the corresponding wind 
properties differ from models without clumping, and check 
whether they are consistent with those derived from our ob- 
servational diagnostics. Remember that if assumption (b) or (c) 
were true, the presently predicted line acceleration is much too 
large. It has to be clarified whether strong clumping is able to 
induce such a large shift in the ionization balance (see Sect. l3.4> 
that the bulk of the accelerating lines are shifted away from the 
flux maximum, such that a reduction in the acceleration is pos- 
sible. 

Finally, time-dependent hydrodynamic simulations must 
also continue. In particular, differences between self-excited 
and triggered structure formation have to be investigated, and 
conditions found which might allow for a much more strongly 
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clumped radio domain than presently predicted (implied if as- 
sumptions (b) or (c) were true). 

In this context, the following, concluding remark is rele- 
vant. Though the usual interpretation of clumping relies on a 
relation to the intrinsic instability of radiative line-driving, the 
issue of whether the redistribution of wind material occurs pre- 
dominantly on small (~0.01 R*) or large (~1 R*) spatial scales 
has not yet been resolved. Small-scale clumping is suggested 
by observations of emissio n-line micro-v ariability in one of 
our targets (HD 66811; see Eversberg et al...l998.) . However, 
structuring of hot-star winds on large scales is indicated by 
the ubiquitous presence of recurrent wind profile variability in 
the form of discrete abs orption component s (DACs; see, e.g., 
IPrin ia & Howarth 19861 iKaper et aljfl996l) . Since there is no 
consensus on the physical origin of DACs, the structure re- 
sponsible for them is not included in the present generation of 
models. 

Future studies will help to address this issue by determin- 
ing whether objects with particularly well-studied DACs (e.g., 
HD 24912, HD 203064, HD 210839) can be modeled success- 
fully without including large-scale structure. The presence of 
unexplained residuals from our self-consistent models (which 
cannot be discounted, due to missing far-IR observations, and 
which might akeady have been identified in the mid-IR fluxes 
of HD 24912, or in the somewhat discordant mm-observations 
of HD 210839, cf. Appendix IbI with small-scale clumping 
would imply that large-scale structures also play a role in the 
redistribution of wind material, and would help to address the 
issue of whether DACs represent localized enhancements in the 
mass flux. 
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Appendix A: The journal of the VLA observations 

is given in Table lA. fi fsee Sect. l2.3> . 

Appendix B: Comments on individual objects 

In the following, we will give, where necessary, some comments on 
the fits for the invidual objects. All results have been summarized 
in Table Q The fits for objects with H„ in emission are displayed in 
Rigs. lB.ll and lB.2l for objects with "intermediate" H^, profile types in 
Rig. lB.3l and for objects with H„ in absorption in Rigs. lB.4l and lB.5l 

B.1. Objects with Ha in emission 

Cyg OB2#7. Ror the hottest object in our sample, only upper limits 
for the radio fluxes are available. The derived mass-loss rate is conse- 
quently an upper limit as well (and the clumping factors corresponding 
lower limits), and based on the assumption that this star is a thermal 
emitter. By means of our regression (Eq. |4}, helium is predicted to 
remain doubly ionized throughout the entire wind (this is the only ob- 
ject in our sample for which this is so), whereas specific models within 
our grid (located in the relevant parameter range) indicate that helium 
might still recombine in the outermost, radio-emitting region. Thus we 
have derived two solutions for this object, both for an ionized and a 
recombined radio regime." 

Ror the doubly ionized solution, we derive a (maximum) mass- 
loss rate of 2.8-10"* Mg/yr. The lower wind is strongly clumped to a 
similar degree in regions 2 and 3 (f^" = 10 and f^"^ = 8. . .12, respec- 
tively. The lower value for /""' results in a good fit of the 10 yum flux, 

" The IR fluxes have been synthesized with doubly ionized helium 
in both cases, since they form well below the radio photosphere. 
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Table A.l. Journal of the VLA observations, including obser- 
vation dates, observing frequencies, time on targets, calibrators 
for flux-density bootstrapping and VLA configuration. 



Star 


date 


freq. 
(GHz) 


time 
(min) 


cal 


conf. 


CygOB2#7 


Feb 15, 2004 


4.86 


60 


3C286,3C48 


CnB 


CygOB2#7 


Feb 15, 2004 


8.46 


45 


3C286,3C48 


CnB 


CygOB2#10 


Feb 15, 2004 


4.86 


20 


3C286,3C48 


CnB 


CygOB2#10 


Feb 20, 2004 


4.86 


40 


3C286,3C48 


C 


CygOB2#10 


Feb 15, 2004 


8.46 


20 


3C286,3C48 


CnB 


CygOB2#10 


Feb 20, 2004 


8.46 


25 


3C286,3C48 


C 


CygOB2#10 


Feb 15, 2004 


14.94 


20 


3C286,3C48 


CnB 


CygOB2#10 


Feb 20, 2004 


14.94 


40 


3C286,3C48 


C 


CygOB2#ll 


Feb 15, 2004 


4.86 


30 


3C286,3C48 


CnB 


CygOB2#ll 


Feb 15, 2004 


8.46 


30 


3C286,3C48 


CnB 


CygOB2#ll 


Feb 15, 2004 


14.94 


30 


3C286,3C48 


CnB 


HD 14947 


Apr 04, 2004 


4.86 


30 


3C48 


C 


HD 14947 


Apr 04, 2004 


8.46 


30 


3C48 


C 


HD 14947 


Apr 04, 2004 


14.94 


30 


3C48 


C 


HD24912 


Mar 09, 2004 


4.86 


20 


3C286 


C 


HD24912 


Mar 09, 2004 


8.46 


15 


3C286 


C 


HD 24912 


Mar 09, 2004 


14.94 


20 


3C286 


C 


HD24912 


Mar 09, 2004 


43.34 


20 


3C286 


C 


HD 34656 


Apr 04, 2004 


4.86 


40 


3C48 


C 


HD 34656 


Feb 09, 2004 


8.46 


45 


3C147,3C48 


CnB 


HD 34656 


Mar 09, 2004 


14.94 


20 


3C286 


C 


HD 34656 


Apr 04, 2004 


14.94 


20 


3C147,3C48 


C 


HD 36861 


Feb 09, 2004 


8.46 


40 


3C147,3C48 


CnB 


HD 36861 


Feb 09, 2004 


4.86 


40 


3C147,3C48 


CnB 


HD 37043 


Mar 09, 2004 


4.86 


30 


3C286 


C 


HD 37043 


Mar 09, 2004 


8.46 


30 


3C286 


C 


HD 37043 


Mar 09, 2004 


14.94 


30 


3C286 


C 


HD 190429 A 


Mar 01, 2004 


4.86 


20 


3C48 


C 


HD 190429 A 


Mar 01, 2004 


8.46 


20 


3C48 


C 


HD 190429 A 


Mar 01, 2004 


14.94 


20 


3C48 


C 


HD 190429 A 


Feb 26, 2004 


43.34 


20 


3C48 


C 


HD 203064 


Mar 01, 2004 


4.86 


60 


3C48 


C 


HD 203064 


Mar 01, 2004 


8.46 


60 


3C48 


C 


HD 203064 


Apr 04, 2004 


14.94 


40 


3C48 


C 


HD 207198 


Feb 20, 2004 


4.86 


60 


3C286,3C48 


C 


HD 207198 


Feb 20, 2004 


8.46 


60 


3C286,3C48 


C 


HD 207198 


Feb 20, 2004 


14.94 


60 


3C286,3C48 


C 


HD 209975 


Feb 20, 2004 


4.86 


30 


3C286,3C48 


C 


HD 209975 


Feb 20, 2004 


8.46 


30 


3C286,3C48 


C 


HD 209975 


Feb 20, 2004 


14.94 


30 


3C286,3C48 


C 


HD 210839 


Feb 26, 2004 


4.86 


20 


3C48 


C 


HD 210839 


Feb 26, 2004 


8.46 


20 


3C48 


C 


HD 210839 


Feb 26, 2004 


14.94 


20 


3C48 


C 


HD 210839 


Feb 26, 2004 


43.34 


20 


3C48 


C 



but slightly too narrow wings of H„, whereas with /""' = 12 we can 
fit these wings perfectly, but so mewhat overes timate the 10 yum flux. 
As for the unclumped models JMokiem et al.ll2005.V the absorption 
trough cannot be fitted well by models with p < 0.9 (nebular emis- 
sion?), though the wings are nicely matched. If we assume, on the 
other hand, that the trough is refilled by the wind alone, the complete 
profile can be reproduced with y6 ~ 1 and f^ = 8, Z™** = 10. . .12, re- 
spectively. From the shape of the trough we derive Hn ~ 11, otherwise 
it becomes too narrow or too deep. 

The alternative solution with helium recombined in the radio re- 
gion yields a considerable larger mass-loss rate, M= 4-10"*MQ/yr, 



since we have adopted a large helium content, i'Hi;=0.21 (compare 
with the case of ^ Pup; see Sect. l4.lt . All clumping properties scale 
accordingly, and the best solution (foryS=0.9) is obtained with f™ = 5 
and f^"^ = 4. . .6. Since the 10 yum flux indicates that helium is not 
completely ionized, even in the outermost IR photosphere (otherwise 
it would lie somewhat higher), we prefer the recombined model for 
our final solution (see Table 0. In the corresponding fit diagram, we 
have indicated both possibilities though (solid: recombined; dotted: 
ionized). 



HD 190429A. For this object, there are two measurements at 3.5 
cm which are co nsiderably different, n amely 200 fi] (our observations) 
and 280 ^J from lScuderi et al] il998h . As is obvious from the fit di- 
agram, the 6 cm flux (our measurement) is consistent with the un- 
published 3.5 cm value provided by Scuderi et al., whereas it lies too 
high with respect to our 3.5 cm measurements. Thus, either the star is 
strongly variable, or a non-thermal emitter, or the errors estimated for 
our observations are too optimistic. Note that the 0.7 cm measurement 
(upper limit) is consistent with our 3.5 cm flux. A "wrong" assumption 
concerning the He recombination cannot explain this dilemma: if the 
ionization degree was higher than predicted, the 0.7 cm flux would be 
most affected and would lie at a level higher than actually observed. 

On the assumption that we see thermal emission and that the dis- 
crepancy is due to measurement problems, the maximum mass-loss 
rate is constrained to lie between 7.5 (dotted) and 9.5-10"*MQ/yr 
(solid), and both limits have been indicated in Table0 By adjustment 
of the clumping factors, we obtain a perfect fit for H„. If the 0.7 cm 
flux is not much lower than its upper limit, f°"^ must be lower than, or 
equal to, 2. The only other discrepancy found for this object concerns 
the 4.63 yum measurement from sCastor & Simon 1 1983) . which cannot 
be matched by any of our models. 



HD 15570 can be fitted without any problems, and the only com- 
plication arises because of the large error bars attributed to the scuba 
fluxes. Since, for the corresponding wavelength, it is not completely 
clear whether He is already recombined or not, we have investigated 
both possibilities. In the recombined case (which is consistent with our 
predictions: recombination at 6.3 R„ 1.3mm radiation becoming opti- 
cally thick at 9.5 R,), the wind must be significantly clumped in region 
4 (/™' = 5. . .20); larger values can be excluded from the H^ wings. If, 



on the other hand, the wind is not recombined in the 1.3 mm-forming 
region, a value of f°^' = 1 is still consistent with the limit of the scuba 
data. The clumping in region 3 (/""' = 6) is well-constrained from the 
Ho, line wings, though a lower value, Z™"" = 4, results when we force 
the 10 yum flux to be matched. In the latter case then, H^, becomes a bit 
too narrow. 

Note that the two measurements at 3.5 cm almost overlap (but 
not completely, indicating a certain variability), and we have forced 
our solution to comply with their average value. In the fit diagram, 
we have plotted three solutions which are consistent with the error 
bars for the scuba measurements: /™' = 5 (dotted), /™' = 13 (solid) 
and f°^" = 20 (dashed). To find even closer constraints on the outer 
wind clumping requires lower error bars. Additional far-IR observa- 
tions (though being important as consistency checks) will not help to 
improve this uncertainty, since the far-IR is insensitive to any reason- 
able variation of /™' for this object. 



HD 6681 1 has been already discussed in some detail; see Sect. 14. II 
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Fig. B.l. Fit diagrams (left: Hq, profile; right: IR/radio continua) for objects with H(j, in emission. Arrows indicate upper limits. For 
parameters, see Table0 Alternative solutions (dotted, dashed) are discussed in the comments on individual objects in ApoendixlBl 
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HD 1 4947. The 3.5 cm flux is well determined (with some variabil- 
ity), whereas only upper limits are available at 2 and 6 cm. The result- 
ing mass-loss rate is M = 8. . . 12 -10"* Mg/yr, and H„ can be perfectly 
fitted, with rather low clumping factors in the lower wind. In the fit di- 
agram and Table we have indicated the intermediate solution with 
M = 10-10-'^ Mo/yr and clumping factors /j," = 3.1 and f^"^ = 2.5. 

Cyg 0B2#1 1 has similar clumping properties to HD 14947, and the 
maximum mass-loss rate can be derived to within small errors: M = 
5 ± 0.5 -10"^ Mo/yr. From H„, the potentially unclumped region must 
be located within ri„ ^ 1.2. From the line wings, Z™"* is somewhat 
larger than /J", and f°"^ might be tightly constrained if far-IR obser- 
vations were available. Problems for this object concern the blue side 
of th e Hg emission being predicted as too narrow, and the 10.9 pm 
flux JLeitherer et alJil982') . which cannot be matched by any of our 
models. 



HD 21 0839. Though the error bars for the scuba fluxes are smaller 
than for HD 15570, two different, barely overlapping fluxes have been 
measured, which might introduce a twofold solution for region 4, 
though the maximum mass-loss rate is well defined. 

A first solution (solid) can be derived for the lower scuba mea- 
surement, with constant clumping in the lower wind, /^!" = /""' = 6.5 
until rout =10, and no clumping in the outer part. The upper scuba 
measurement can be fitted by additional clumping in region 4, with 
/™' = 5. . . 20, but in this case the 0.7 cm flux appears as too large. 

A second, slightly better solution (which is indicated in Table Q 
can be found if one assumes constant clumping (again with /^ = 6.5) 
until r < 4, and a larger clumping factor of fa = 10 until r < 15. 
With /™' = 1, the lower 1.3 mm flux is matched (dotted), whereas 
with /™' = 8 the upper one can be fitted (dashed). As before, however, 
the 0.7 cm flux is then predicted as too large. For all solutions, H„ is 
perfectly reproduced, and a value of r,„ S 1.2 can be constrained from 
its trough. 

HD 1 92639. Only one radio measurement is available, and only as 
an upper limit. Adopting this value and assuming thermal emission, 
the maximum mass-loss rate can be restricted to M ^ 3- 10"* Mo/yr, 
with constant clumping factors, /^™ = /""' = 3.5, in the lower wind, 
and Hn ^ II- F or all our simulations, the observed 4.63 /im flux (taken 
from lCastor & Simon 1983) is smaller than synthesized, though better 
reproduced than for HD 190420 A, and independent of the ionization 
equilibrium for helium. 

HD 3061 4 is perfectly matched, both in the radio and in H„, with a 
moderate degree of clumping in the inner and intermediate wind. 



this measurement. Only low clumping factors are required to fit H^,, 
though higher values would be necessary if the mass-loss rate were 
lower. Note that with f^^"^ = 2.0 the wings of H„ are nicely matched, 
but the 10 yum flux is slightly overestimated. With /J™'' = 1, on the 
other hand, the latter problem can be cured, at the expense of H„. 



Cyg 0B2#1 can be fitted accounting for weak clumping in the 
lower wind (/j" = 1.4, /J™'' = 1.8) if j8 is left at its nominal value 
of 1.05. Hn must be S 1.2, and clumping effects are seen only in the 
inner wind. The observed 10 yum flux is larger than predicted, which 
cannot be corrected for by a non-recombined wind, as the temperature 
is too low for such a scenario. 



6.3. Objects with Ha in absorption 

For all objects with H„ in absorption, we have used rout = 10, since 
due to the lower wind density, the IR and radio emission is formed at 
smaller distances from the star (cf. Sect. 14. H . E.g., for HD 36861, the 
wind becomes optically thick at 2 cm only for r S lORi,. 



Cyg OB2#8C remains rather unconstrained by our analysis, since 
only one upper limit in the radio range is available (at 6 cm), and this 
upper limit yields a mass-loss rate larger than the one derived from H„ 
(M= 3.5- 10"* Mo/yr for/3 = 1). Thus, the largest possible mass-loss 
rate has been adopted from this value, and the only definite statement 
concerns f™'' being similar to f^". 

HD 34656 is the only object within our sample where the radio 
mass-loss rate (if thermal emission) is definitely larger than the H^ 
mass-loss rate.^" Unfortunately, only one measurement (at 3.5 cm) 
provides a hard number, whereas the 2 and 6 cm measurements yield 
upper limits only. Thus, non-thermal emission cannot be excluded, 
at least to some extent, if one compares the 3.5 and 6 cm fluxes 
(Fig. IB.41 . Besides being a non-thermal emitter, there are two other 
possibilities: either the 3.5 cm measurements are somewhat corrupted 
(i.e., can be regarded as upper limits only), or the outer wind is more 
heavily clumped than the inner one. In the latter case, the maximum 
mass-loss rate results from H„ instead of from the radio, and corre- 
sponds to 3-10~* Mo/yr for/3 = I, which is the lowest possible value 
such that /? remains consistent with our data (wings of H„). Note that 
a value of /3 = 1.1 and M = 2.6- 10"* Mo/yr, as derived in Paper I, 
gives a slightly better solution. Since the (thermal) radio mass-loss 
rate corresponds to a value of 7-10"* Mo/yr (and would result in an 
emission profile for H„; cf. the dotted solution), clumping factors of 
yftr _ (M,,,;jio/MHff)" ~ 6 are necessary to obtain a simultaneous fit, 
with /^■j'"' = /^'j" being well-constrained. 



B.2. Objects witli "intermediate" Ha profile type 

Cyg OB2#8A js a confirmed non-thermal radio emitter 

JBiesing et alJ 1989^. In order to obtain at least an estimate, as 
low a maximum mass-loss rate as possible has been adopted (from 
the 2 cm flux), although this might still be even smaller, of course. 
With jS = 0.74 (taken from the optical analysis using homogeneous 
models. iMokiem et al. 2005 ), the wings of H„ are fitted best, whereas 
the absorption becomes too deep. A value of j8 = 0.85 (dotted) 
improves the trough, but the emission then becomes too large. The 
20 yum flux indicates that our prediction for the recombination radius 
of helium might be erroneous, and a completely recombined model 
(which at these temperatures is rather improbable) can indeed fit 



HD 2491 2. All radio fluxes measured for^ Per are upper limits. Our 
fit diagram shows that the 0.7, 2 and 6 cm limits, if taken at face 
value, are consistent with thermal emission and a maximum mass- 
lo ss rate of 2.3- 10~* M^/ yr, which is very close to the value provided 
bv lRepolust et alJi2004) . using unclumped models. When accounting 
additionally for the 3.5 cm flux, one derives a maximum mass-loss rate 
of 1.2- 10"* Mo/yr. In the following we will consider both possibilities. 
The solution with larger M requires weak clumping in the lower- 
most wind (/j!" = 2.1), and additional clumping in region 3 (/^'J'"' = 
5), if the small emission humps on the red and blue side of the H„ ab- 
sorption are due to clumping and not to other processes (see below). 



For HD 209975, both mass-loss rates overlap within the errors. 
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Fig.B.3. As Fig. IB. II but for objects with "intermediate" H„-profile type. 



The maximum value of /™' is restricted by f°"' i 2, otherwise the 
radio band becomes affected and the maximum mass-loss rate must 
be decreased. 



HD 203064. For the standard value of /3 = 0.9, a model with 
all clumping factors being unity is consistent with the observations. 
Values for /'!"'' S 2 can be excluded. 



The lower M solution (which is consistent with all radio measure- 
ments) requires considerable clumping in the lower wind. Assuming 
a "standard value" of f} = 0.9 (which has been used for most of the 
following objects as well, but see Sect.|43j, f^ = 8 and f^"^ < 20..25, 
the humps can be explained by clumping. Furthermore, the unclumped 
region (if any) can be constrained by ri„ S 1.1. For this solution, 
/™' S 3, otherwise M is even lower. In our fit diagrams, we have 
plotted the high M solution (solid), the low M solution (dotted) and 
the low M solution with f?"' = 1, which does not fit the emission 
humps in H„. At least this uncertainty might be resolved if future far- 
IR measurements become available. 



Let us comment finally on the strong excess measured in the mid- 
IR, between 8.7 to 11.4 jixa (taken from Gehrz et al. 197 4), which is 
in stark contrast to the 12 //m IRAS data from Beichman e t al. ( 1988). 
This discrepancy (see also Sect.|6j cannot be due to a wrong flux cali- 
bration, since measurements from the same source have been used also 
for HD 30614 and HD 36861, without any apparent problems. Thus, 
^ Per is either strongly variable in the mid-IR, or the mid-IR excess 
is due to smother physic al process (e.g., co-rotating interaction zones, 
see lde Jong et all200lL or a wind compressed equatorial region). The 
latter interpretation in particular is consistent with the red and blue 
emiss ion h umps observed in H„, which have also been seen in Hen 
4686 iHerreroetaljl992. Fig. 4). 



HD 36861 . All radio measurements provide only upper limits, and 
we have indicated a model with a consistent maximum mass-loss rate, 
M = 0.4-10"* Mg/yr. For this value and/3 = 0.9, the innermost clump- 
ing is weak again, /'," = 2 



lo/yr. For this value and/3 = 

. 4, and /""' must be lower than 20 (from 
the wings of H„). Solutions with /"j"' 



> 2 are no longer consistent 
with the adopted mass-loss rate. In Fig. IB. 51 we have indicated the 

1, solid) and maximum 
(/^»i = 4, ff'^ = 20, dotted) clumping. Note that only far-IR or mm 
observations will help to disentangle this uncertainty. 



solutions with minimum (/^!" = 2,/^ 



mid 
cl 



HD 2071 98 has well-defined radio measurements, and an un- 
clumped wind with M = 1.0. . . 1.2-10"'' Mo/yr (for^S = 0.9) matches 
all observational constraints (f'™'^ S 2). The 2 cm flux can be re- 
produced with f°^' = 10 (and <15); larger values are excluded by 
the H„ wings. Displayed are the solutions for an unclumped wind at 
M = 1.0-10"* Mo/yr (solid), a wind with additional clumping in the 
outer region (f°^' =10, dotted) and a homogeneous wind with M = 
1.2-10"" Mo/yr (dashed). 

HD 37043. By inspection of the measured radio fluxes, this star is 
either a non-thermal emitter (SB2!), or the 6 cm flux is erroneous. At 
3.5 cm, we have two measurements which are consistent. To obtain 
more conclusive results, one needs to re-observe this star in the radio 
range. 



1.05 r^ 



J. Puis et al.: Constraints on the clumping factors in hot star winds 

CygOB2#8C CygOB2#8C 



37 




6500 6520 6540 6560 65B0 6600 6620 6640 
wavelength (A) 




10 



HD34656 



10^ lO'- 10^ 

wavelength [/xm] 

HD34656 



10^ 10" 



1.2f 




OJ 0,9 r 



O.H h 

6500 6520 6540 6560 6580 6600 6620 
wavelength (A) 




10 



HD24912 



10" lO'' 10" 

w^avelength [urn] 

HD24912 



10" 



h05 



1.00 '", 



0.95 - 



0,90 



0.85 - 




1.00 -i^ 



?f 0.90 - 



0.80 

6500 6520 6540 6560 6580 6600 6620 6640 
w^avelength (A) 



10" 10'' 10" 

wavelength [,uin] 



Fig.B.4. As Fig. IB. II but for objects with Hq, in absorption. 
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Fig.B.5. AsFig.lMl 



Nevertheless, we present two solutions: an upper one discarding 
the 3.5 cm data and being consistent with the upper limit at 2 cm, and 
a lower, more likely one (which would also be an upper limit if the 
object were a non-thermal emitter), discarding the 6 cm measurement. 

In the first case, M = 0.8- IQ-* Mo/yr, and a smooth wind is consis- 
tent with the observations. In the second case, M = 0.25-I0-* Mo/yr, 
and the wind is strongly clumped at least in the lowermost wind, with 
/^'" = 12 for Hn = 1. . .1.05 to r^id S 1.3. Due to the low density, 
clumping in other regions has a very low impact on the model fluxes, 
and we can exclude only values /""* > 20 and /^™" > 10 (otherwise 
the maximum mass-loss rate must be lower). Plotted are the "smooth" 
solution with the upper value for M (solid), and the lower M solution, 
which is strongly clumped (dotted). 



HD 209975 has already been discussed in Sect. l4.1l 



